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Abstract 

We prove the asymptotic stability and asymptotic completeness of an arbitrary number of weakly 
interacting solitons for NLS. 

1 Introduction 

The nonlinear Schrodinger equation 

(NLS) i^- = -Aip - (3(\tp\ 2 )ip, xeW 1 

has in general (exponentially) localized solutions in space, provided the nonlinearity has a negative 
(attractive) part. This is due to a remarkable cancellation of the dispersive effect of the linear part with 
the focusing caused by the attractive nonlinearity. To find such solutions, we look for time periodic 
solutions ip = e %ult <j> u {x). It follows that (f> u , if it exists, is a nonzero solution of the problem 



(ENLS) -u<f> u = -A<f> u - /3(|^| 2 )^. 



We shall refer to such solutions as nonlinear eigenf unctions. In general, for <f> u to be localized (at least 
as I? function) we need u> > 0. 

The general existence theory for this elliptic problem has been studied in great detail, see the work 
of Coffman [Cof], Strauss [Str], and Berestycki, Lions [BL]. 

It is easy to see that if ip^itjx) = e lujt <j3 UJ (x) is a solution of NLS, then for any vector a £ R n the 
function (j>u(t,x — a) is also a solution. More generally, NLS is invariant under Galilean transformations 

(1.1) ftr.uCt) := e -^-5M 2 y(^+ B )P, 

and therefore we can construct solutions from <j> u which are moving with arbitrary velocity v. As a 
result we obtain a family of exponentially localized solutions 

^, 7 ,A" = e^-^-^^Ux -vt-D) 
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parametrized by a constant (2n + 2) -dimensional vector (v,j,D,uj), which are known as solitons. 

Solving the initial value problem for NLS requires understanding of two fundamental questions. 
The first is the existence of global in time solutions. Due to the focusing character of the nonlinear 
term the global existence theory ought to be based on the 1? conservation law. By the results of 
Kato [Kato] and Tsutsumi [Tsu] we can construct unique global solutions for any I? initial data under 
the assumption that the nonlinearity (3 satisfies the condition \/3(s)\ < (1 + \s\) q with q < ^ . 

The second problem is that of the asymptotic behavior of solutions as t — > +oo. Guided by the 
completely integrable models in one dimension we expect to have solutions with the asymptotic profile 
of N independently moving solitons: 

N 

(1-2) l/>(t,x) ~^2^v k , lk ,D k ,uj k (t,x). 

k=l 

Moreover, given such a solution ip(t,x) we expect that initial data ip(0,x) + R for a suitably small 
perturbation R leads to a solution ip with \\il>(t, •) — ip(t, -)|| — > as t — > oo in an appropriate norm. 
The latter property is known as asymptotic stability. In the context of NLS one needs to modify (1.2) 
since one needs to make (v k , j k , D k ,uj k ) time-dependent. This fact was already observed in the context 
of orbital stability, see Weinstein [Wei]. 

In this paper we give an affirmative answer to the question of existence and asymptotic stability of 
solutions with iV-soliton profiles under the assumption of weak mutual interaction between the solitons. 
A superposition of TV spatially separated moving solitons gives only an approximate solution of NLS. 
Our goal is to show that the initial data 

N 

(1-3) Mx) = E e Wfc(0) "^ (0) &, fc (o)(s " D k (0)) + Ro, 

k=l 

give rise to an asymptotically stable solution with the profile of N independent solitons, with perhaps 
different parameters (4,7fc, D k ,uj k ) / (4(0), 7fc(0), D k (Q), uj k (0)). 

The function Rq is a perturbation satisfying a smallness assumption on its L 1 n L 2 norm together 
with its derivatives. An important aspect of our main result is the assumption that the solitons are 
weakly interacting. This condition can be enforced in two ways. Firstly, one can assume that the 
initial shifts D k (0) and the initial velocities v k are chosen to model the case of non-colliding solitons: 

(1.4) \D k + v k t - D k < - v k d\ > L + ct, k^k' 
for some sufficiently large constant L. 

Alternatively, one can assume that the relative initial velocities of the solitons are large, i.e., 

(1.5) min \vj — v k \ > L. 

For the most part, we give details only for the case of (1.4) and leave the simple modifications required 
by (1.5) to the reader. Note that (1.5) does not rule out that the solitons collide. However, in view 
of (1.5) the time of interaction is of size L" 1 , and therefore the overall interaction remains weak. 
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We shall also require certain spectral assumptions on the nonlinear eigenstates 0^ which will be 
explained below. 

We believe that the methods we use may be applied to other classes of equations with solitary type 
solutions and other symmetry groups (e.g. Lorentz instead of Galilean), if and whenever certain linear 
L p decay estimates can be verified for the linearized operators around one such soliton. A detailed 
analysis of such L p estimates for NLS was recently given in [RSS]. 

To explain our results we recall the precise notions of stability. Suppose we take the initial data of 
NLS to be an exact nonlinear eigenstate <j) w , plus a small perturbation Rq. What is then the expected 
behavior of the solution? If the solution ip(t) stays near the soliton ^(t) = e lu ' t (f) u ,(x) up to a phase 
and translation for all times (in H 1 norm) we say that the soliton ip u is orbitally stable. If, as time 
goes to infinity, the solution in fact converges in L 2 to a nearby soliton plus radiation 1 we say that the 
solution is asymptotically stable. 

Orbital stability of one soliton solutions has been subject of extensive work in the last 20 years. 
The first results date back to the work of Cazenave and Cazenave-Lions on logarithmic and monomial 
nonlinearities. The general case has been treated in the defining works of Shatah-Strauss [ShSt], 
Weinstein [Wei] and [We2], and Grillakis-Shatah-Strauss [GSS1]. The general phenomena that has 
emerged from their results is that the orbital stability is essentially controlled by the sign of the 
quantity dull^llz, 2 (stable, if positive, and unstable, if negative). However, all these results addressed 
orbital stability of a special class of solitons generated by ground states: positive, radial solutions of 
the equations ENLS. 

In [BL] Berestycki-Lions proved the existence of a ground state in three or more dimensions for any 

9 4 

to ^ under the conditions that the nonlinearity (3 verifies lim (3(s )s n ~ 2 > and such that there 

exists < so < oo, with G(sq) > 0, for G(s) = 2J^((3(t 2 )t — lut)cIt. In fact, in their work ground 
states are found as minimizers of the constrained variational problem: 

(1.6) inf J[u] =inf ( f |Vu| 2 : f G{u) = l) . 

The question of uniqueness of a ground state has been studied in [McS], [Kw], [McL]. 

The asymptotic stability of one ground state soliton solutions of NLS and other equations was first 
shown for NLS with an extra attractive potential term in [SW1], [SW2], and [PW], for one NLS soliton 
in dimension one in [BP1] and in dimensions n > 3 in [Cu]; for NLS-Hartree see [FTY]. 

While the arguments for orbital stability were essentially based on Lyapunov type analysis and 
relied only on some limited information about the spectrum of the associated linear problem, the 
proofs of asymptotic stability required much more detailed properties of the related linearized systems. 
In particular, it led to the need to impose additional spectral assumptions on the linear operators 
associated with a soliton. 

To describe the problem of linear stability (and spectral theory) we linearize NLS around a soliton 

: a function with asymptotic behavior e ttA f. Alternatively, we can and will replace the L 2 by the L°° convergence. In 
the latter topology the contribution of the radiation can be ignored. 
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w = e J Vo>> using the ansatz 
(1.7) ip = e w [(f) U) (x-vt) + R(t,x-vt) 



The resulting "linear" operator acting on R, also has a term containing R. After complexifying the 
space to (R, R) we are left with a matrix non self-adjoint operator of the type 



(1.8) H 
acting on I? x I? . Here 



L+ W(x) 
-W{x) -L+ 



L+ = _ A + u - (3(4>l) - p {<&)<&, W(x) = -/?(<&)<& 

The operator L + is a self-adjoint perturbation of —A by an exponentially localized function, and W 
is an exponentially localized potential. Moreover, centering the perturbation R around the soliton as 
in (1.7) ensures that 7i is time-independent. Note that by non self-adjointness of H one can no longer 
guarantee that sup t ||e* tw ||2->2 < oo where U(t) = e ltn is given by id t U + HU = 0, U(0) = Id. In fact, 
the operator TL has a zero root space M := (J^>i ker?^ of dimension at least 2n + 2 containing the 
eigenfunction (j)^ as well as the elements generated from <j> w by infinitesimal symmetries of the problem. 
We decompose 

L 2 x L 2 =J\f + N* ± 

and let P denote the projector on the second term in this decomposition (here M* = U^>i ker(7i*)^). 
It is easy to see that ||f7(t)/||2 grows polynomially for some / G J\f. On the other hand, it is known 
from work of Weinstein [Wei] that 

(LS) sup || U{t)Pip \\ L 2 < oo 

t 

under certain conditions on the nonlinearity and provided cf)^ is the (positive) ground state of ENLS. 
Generally speaking, we refer to the property (LS) as linear stability. 

In the case when H. has one negative eigenvalue and (p^ is the unique ground state, one can show 
that the condition cr(Ti) c R is equivalent to the orbital stability condition 9a) 1 1 ^ 1 1 l 2 > ( see [Gr], 
[BP1]). Although due to the lack of self-adjointness this is not sufficient for linear stability, additional 
arguments show that (LS) in fact holds just under the above conditions (see [Wei], [GSS1]). 

Linear stability of ground states has been considered for a large class of NLS in the work of 
Weinstein [Wei], [We2], and Shatah, Strauss, Grillakis , [ShSt], [GSS1], [GSS2], [Gr], see also [SuSu], 
[Stu]. However, unconditional results were established only in the case of monomial nonlinearities 
P(s) = s p . Moreover, linear stability has been shown to be essentially equivalent to the orbital stability 
(except in the case of an I? critical nonlinearity (5(s) = s 2 / n ). 

Linear stability plays an essential role in the results on asymptotic stability. Moreover, the proofs of 
asymptotic stability of one-soliton solutions required even more stringent assumptions on the structure 
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of the spectrum of TL, see Buslaev, Perelman [BP1], and Cuccagna [Cu]. This can be linked to the fact 
that on the linearized level asymptotic stability requires dispersive estimates of the type 



(1-9) \Mt)P4o\\ L oo<t- 



^OllLi 



[Cu] or similar L 2 -weighted decay estimates. To prove such estimates one needs to impose additional 
spectral conditions such as: absence of the discrete spectrum for TL on the subspace J\f ± , absence 
of embedded eigenvalues, and absence of resonances at the edges of the continuous spectrum. The 
dispersive estimates for such matrix Hamiltonians in dimensions n > 3 were proved by Cuccagna by 
an extension of the method introduced by Yajima in the scalar case [Yal]. In that approach, the 
decay estimates follows as a consequence of the proof of the L p — > L p , Vp € [1, oo] boundedness of the 
wave operators. In our recent work [RSS] we suggested a perhaps more straightforward approach for 
proving such estimates which instead relies on construction of the analytic extension of the resolvent 
of TL and goes back to the work of Rauch [Rau] in the scalar case. This method, however, requires that 
all potential terms in the Hamiltonian TL are exponentially localized functions, which perfectly fits the 
problem at hand. We refer to Hamiltonians verifying the required spectral assumptions (as well as the 
linear stability condition) as admissible. The proof of the dispersive estimates crucially relies on the 
time independence of the Hamiltonian H, which was ensured by the choice of the ansatz. 

In our work we choose the initial data of the form (1.3), with cf) u}k verifying the elliptic prob- 
lem ENLS, satisfying the separation condition (1.4) and study the time asymptotic behavior of the 
corresponding solution of the time-dependent NLS. 

As in the study of asymptotic stability of one-soliton solutions the first objective is the analysis of 
the linearized problem. It is natural to use the following ansatz for the solution: 

N 

(1.10) il>(t,x) = w + R(t,x) = ^2w k + R{t,x), w k (t,x):=e iek c/) <JJk (x-v k t-D k ) 

k=i 

where 9 k are the phases associated with k-th soliton and v k ,D k ,uJk are its parameters 2 . We substitute 
this ansatz into the equation (NLS) and retain only those terms which are linear in R, R. The resulting 
linear problem for the unknown (R, R) is 

(1.11) id t U + H{t)U = 0. 
It contains a time- dependent complex Hamiltonian 

wm - ( L+ w ^ x < 

n{ ) ~ \-W(t,x) -L+ 
where 

L+ = -A - (3(\w\ 2 ) - p'(\w\ 2 )\w\ 2 , W(t,x) = -(5{\w\ 2 )w 2 . 

Because of the smallness assumption on the initial perturbation Rq we can assume the the pa- 
rameters of the final asymptotic profile will lie in a small neighborhood of the initial parameters 
2 which in the true ansatz are to be made time-dependent 
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(v k {0),'-f k (0), D k (0),uj k (0)). The separation condition on the initial parameters and the exponential 
localization of functions 4> Wk guarantees that w is decomposed into a sum of functions of essentially 
disjoint support. We thus can replace the Hamiltonian H(t) with 



N 



n 19 n -h/(t\ = f~ A Q^+Srft^-Vkt-Dk) W k {x-v k t-D k ) 

■ ' " !! " I A) + 2^{-W k (x-v k t-D k ) - Uk (x - v k t - D k ) 



k=i 



where 



U k = -^l k {x))-(5\^l k {x))4>l k {x)f, 
W k = -e™^{cg k {x))cti k {x) 

This Hamiltonian belongs to the class of the so called matrix charge transfer Hamiltonians. Each of 
the Hamiltonians 

wm _f"A 0\ ( U k (x-v k t-D k ) W k {x-v k t-D k )\_ 

Hk[t) ~ I Aj + {-W k (x - v k t - D k ) -U k (x - v k t - D k )) ~ Ho + V ^ x ~ ^ 

represents a linearization around the fcth soliton. Moreover, a specially chosen Galilei transform of the 
type (1.1) maps a solution of the linear problem idtU + H k (t)U = into a solution of the problem 
idtU + H k U = with a time-independent Hamiltonian H k of the form (1.8). Here Hq = diag(— A, A) 
and V k is an exponentially localized complex matrix potential. 

Once again, at the linear level, the heart of the problem of asymptotic stability of N-solitons are the 
dispersive estimates for the soutions of the equation (1.11) with a matrix charge transfer Hamiltonian 
7i(t). Observe that in contrast to the one soliton case, the linearized Hamiltonian is time-dependent. 

Due to the separation condition the problem (1.11) admits "traveling" bound states generated by 
the discrete spectrum of each of the Hamiltonians 7i k . These bound states are formed by eigenfunctions 
and elements of the root space of 7i k boosted by the Galilei transform corresponding to the parameters 
v k ,D k . The paper [RSS] establishes dispersive estimates 

(1-13) \\U{t)M\L*+Lo° < (l + *r%o|kinL2 

for the solutions of the linear time-dependent Schrodinger equation (1.11) with a matrix charge transfer 
Hamiltonian H(t) of the type (1.12) in dimensions d > 3. These estimates hold under the assumption 
that each of the time-independent Hamiltonians H k is admissible (requiring spectral assumptions and 
the linear stability condition) and that the solution U (t)ipo is asymptotically orthogonal to all traveling 
bound states of H k (t). The latter means that for all k = 1, . . . , N 

\\P b (n k ,t)U(t)^ \\ L 2 ->0, as t^+oo 

with Ph(7{ k ,t) denoting the time-dependent projection on the k-ih subspace of traveling bound states, 
i.e., it is the conjugation of the spectral projection of the stationary operator 7i k onto its bound states 
by suitable Galilei transforms. The estimate (1.13) is our main linear estimate. However, to ensure 
that the perturbation R(t,x) in the decomposition (1.10) is asymptotically orthogonal to the subspace 
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of traveling bound states and thus decays in the linear approximation with the rate of t~z, we need 
to allow the soliton parameters a k = (vk, Dk, 7fc, cok) to become time-dependent. This in turn makes it 
necessary to change the form of Wk(t, x), for example, x — Vkt — Dk becomes 

x - Vk{s) ds - D k (t). 
Jo 

The resulting nonlinear problem for R, after the complexification Z = (R, R), takes the form 
(1.14) id t Z + H(t, o(t))Z = a ■ d a w + jSf(Z, w), 

where H(t,a(t)) is a time-dependent Hamiltonian. It is of the matrix charge transfer type provided 
a(t) =const. The term d a w denotes the derivative of the solitary approximation w(t,x;a(t)) = 
J2k=i w k{t, x; cr(t)) with respect to its parameters Ok, the term a denotes the time derivative of the 
soliton parameters cr^, and J\f(Z, w) is a nonlinear term in Z. We introduce the notion of an admissible 
path o~{t) in the space of parameters and a reference Hamiltonian H(t, a) at infinity corresponding to 
the matrix charge transfer Hamiltonian with fixed constants o~k = o~k(t = oo) = (vk, Dk, 7fc, uJk)- At 
each time t the solution Z{t) is required to be orthogonal to the traveling bound states of the charge 
transfer Hamiltonians H(t,a = a(t)) obtained by fixing the parameters a k = (Jk(t) at a given time 
t. This leads to the so called modulation equations for o^, which couple the PDE (1.14) for Z with 
an ODE for the modulation parameters crk(t). To impose the orthogonality condition we first need 
to verify that it is satisfied initially. Using standard arguments, see e.g. [BP1], one can ensure this 
property by modifying the soliton parameters cx(0) slightly in the decomposition of the initial data ifro, 

N 

We later justify the orthogonality condition at any positive time by showing that it is propagated. To 
handle the nonlinear equation (1.14) we introduce the Banach spaces X s and y s of functions of (t,x) 



\f\\ Xs = sup(||^,-)ll^ + (l+i)^EH Vfe /^-)IU 2+ Loc) 

F\\y, = sup^( /V fc ^,-)||L 1 ^ + (l + t)f+ 1 ||V fe F(i,-)llLO- 



(1.15) 



The space X s is designed to control the solution Z{t) itself, while the space y s takes care of the 
nonlinear terms appearing as inhomogeneous terms in the equation for Z. We rewrite the equation 
(1.14) replacing the Hamiltonian H(t,a(t)) with the reference charge transfer Hamiltonian H(t,a). 
The solution operator of the corresponding linear problem 

(1.16) id t U + H(t,a)U = F 

maps 34 - ► %s for any integer s uniformly in a, provided that the solution is asymptotically orthog- 
onal to all traveling bound state of H(t,a). We then show that the nonlinearity F arising from the 
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inhomogeneous terms in the equation Z, written relative to the reference Hamiltonian H(t,a), maps 
X s —> ys f° r an Y s > [§] + !• Given the smallness assumption on the initial data for Z = (R, R) this 
allows one to conclude the desired properties of R. The modulation equations for a(t) are then used 
in turn to control the path cr(t). In particular, we show that there exist final values of the parameters 
a = cr(oo), thus justifying the introduction of the reference Hamiltonian at infinity. The estimates for 
the inhomogeneous problem (1.16) phrased in terms of the mapping between y s and X s are essentially 
the linear dispersive estimates for the time-dependent Schrodinger equation with a matrix charge trans- 
fer Hamiltonian 7i(t) = H(t,a) proved in [RSS]. As it was mentioned before, these estimates only 
require the admissibility of the corresponding individual time-independent matrix Hamiltonians Tik 
(representing the linearization around each individual soliton). The admissibility conditions can be 
somewhat loosely divided into two categories: 

1) Conditions related to linear stability: sup t ||e rfWfc /||2- + 2 < oo 

2) Spectral assumptions on Hk (absence of the embedded eigenvalues, resonances, etc.). 

Remark 1.1. We also require the absence of "spurious" eigenvalues. This means that we assume that 
all of the discrete spectrum of 7ik is generated purely by the nonlinear eigenfunction aJfe , i.e., it 
is described by the generalized eigenspace space of TCk which has the precise dimension (2n + 2) 
(coinciding with the dimension of the parameter space of a^). This is motivated by the requirement 
that the solution Z has to be orthogonal to all traveling bound states of Hk for the dispersive estimates 
to hold, which can only be achieved by the choice of the parameters a^. In general, it is believed that 
the states corresponding to the spurious eigenvalues decay in time but the mechanism of this decay is 
purely nonlinear. We do not pursue this issue here. 

The key ingredient in establishing the linear stability is known to be the monotonicity or convexity 
condition 

(1.17) <9J<Ml 2 >0 or (L;Va,,^)>0. 

The known examples of when the monotonicity condition can be verified are limited to the case of the 
ground states 4> corresponding to the monomial subcritical nonlinearities 

(1.18) I3{s 2 ) = s p -\ Kp<l + -, ± 

n 

and the nonlinearity of the mixed type (see [Sh] ) 

/3( a 2) = s 2 _ s 4 

for the values of uj close to 3/16. In this paper we find a new class of nonlinearities satisfying condition 
(1.17). These nonlinearities lie "near" the subcritical monomials of (1.18) but vanish much faster near 
s = 0. More precisely we consider functions 

(1.19) (3e(s 2 ) = S?- 1 f{s2) 



+ f(s 2 ) 

with a constant 8 > and the function / satisfying the conditions 

(1.20) C\s r+1 < \f(s 2 )\ < C 2 s r+ \ s 2 \f'(s 2 )\ < C 2 \f(s 2 )\, p e (1, 1 + 1), r > -1 
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and prove that given a sufficiently small neighborhood U in the space of parameters u there exists a 
sufficiently small value 6q sucn that for all 6 < 6q and all lo € U the ground state of (3g corresponding 
to a satisfies the monotonicity condition. 

We note that the higher rate of vanishing of (3(s 2 ) at s = is important for asymptotic stability. In 
particular, it should be mentioned that if the power p in the monomial example is too low (p < 1 + ^) 
even the scattering theory (asymptotic stability of a trivial solution) fails. 

We now describe the structure of the paper. 

Section 2. contains the statement of the main result together with the definitions of some of the 
fundamental objects used in the proof. The latter include the definition of an ansatz ip = w a + R, 
separation condition on the initial data ensuring that the solitons are only weakly interacting, the 
notion of an admissible parameter path cr(t), and the spectral assumptions. It also contains the first 
discussion of the linearized Hamiltonians appearing in the later sections. We note here that traditionally 
the results on the asymptotic stability require the smallness assumption on the initial data in weighted 
Sobolev space. Our result uses instead the Sobolev space based on the intersection of L 1 n L 2 (R n ), 
which has a distinct advantage of being translation invariant. 

Section 3. gives a detailed description of the linearization of the equation (1.11) around an 
N-soliton profile w a and introduces the notion of the reference charge transfer Hamiltonian at infinity. 

Section 4. describes the structure of the nullspaces of the Hamiltonians Hj(a), associated with 
the linearization on each individual nonlinear eigenfunction <pj. 

Section 5. recalls the dispersive estimates for solutions of the time-dependent Schrodinger equa- 
tion with a charge transfer Hamiltonian. 

Section 6. derives a system of ODE's for the modulation parameters a{t) = (a±, . . . ,<jn) with 
&k = (vk, Dk,jk,ctk) by requiring the complexified perturbation Z = (R, R) to be orthogonal to the 
the unstable manifold comprised of the elements of the nullspaces of Hj(cr(t)). 

Section 7. gives the bootstrap assumption on the size of the perturbation Z and the admissible 
path a(t) and provides estimates on the difference between the linearized Hamiltonian H(a(t)) and 
the reference Hamiltonian H(a, t) at infinity. 

Section 8. provides the solution of the modulation equations for the path cr(t). 

Section 9. solves the nonlinear equation for the complexified perturbation Z in the space X s for 
s > [^] + 1. This includes algebra estimates designated to show that the JVnorm of the nonlinear 
terms in the equation for Z can be controlled by the ^f s -norm of Z itself. 

Section 10. proves the scattering result showing that as t — > oo the solution ijj(t) decomposes into 
the sum of solitons and a solution of the free linear Schrodinger equation. 

Section 11. discusses the existence for the coupled PDE-ODE system for Z and a(t). 

Section 12. returns to the detailed discussion of the associated linear problems and proves some 
of the assertions made in the first part of the paper. We give the precise definition of an admissible 
Hamiltonian and start the investigation of their spectral properties. In particular, in Section 11.2 we 
describe the spectrum of an admissible Hamiltonian and prove exponential decay of the elements of 
its generalized eigenspaces. In Section 11.3 we specialize to the admissible Hamiltonians arising from 
linearization around a nonlinear eigenfunction </>. We introduce and discuss the associated self-adjoint 
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operators L + and L_ and show that the admissibility conditions on our Hamiltonian (excluding the 
assumptions on absence of the embedded spectrum and resonances) can be reduced to the monotonicity 
condition (1.17), and the statement that the null space of L_ is spanned by <p while L + has a unique 
negative eigenvalue with the corresponding null space spanned by d x d>. Some of the arguments in this 
section follow those of [Wei], [BP1]. 

Section 13. shows that for a very large class of nonlinearities the operator L+, obtained by 
linearizing on a ground state, has a unique negative eigenvalue. 

Section 14. establishes the desired properties of the operators L + and L_ for a particular class 
of nonlinearities (namely those used in our main result) in the case when the nonlinear eigenfunction 
d> is a ground state. 

For the most part of the paper (until Section 12) we do not specify the nature of nonlinear eigen- 
states 4>u k of the elliptic problem (ENLS). In particular, we do not require them to be ground states. 
Instead, we choose to formulate a more general conditional result dependent upon verification of certain 
precise properties of the linearized operators TLk associated with each (j> Uk . It is only in Section 12 that 
we verify some of these assumptions in the case when <fi Uk are ground states. The reason for choosing 
this approach is to emphasize the method which allows us to handle weak interactions of (non-colliding) 
solitons generated by d> UJk , provided that certain properties of each individual eigenstate <p Wk hold. We 
believe that our method will have an even wider range of applications than described here. 

2 Statement of results 

Consider the NLS 

(2.1) ^V> + ^AV + /3(H 2 )V = 
in M. n , n > 3, with initial data 

N 

(2.2) Mx)=Y,™j(Q^)+Mx)- 

Here Wj(0,x) are nonlinear eigenfunctions generating the solitons 

(2.3) Wj (t,x) = w(t,x;a j (0)) = e ie i (t > x) d>(x - Xj(t), ctj(0)) 

(2.4) 9j{t,x) = « J -(0)-x-i(|t; j (0)| 2 -^(0))* + 7 J -(0) 

(2.5) Xj (t) = vj^t + DjiO). 

and <j) = (f>(-, a) is a solution of 

(2.6) |a0-^0 + /5(|0| 2 )0 = O. 
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The solitons Wj as in (2.3) satisfy (2.1) with arbitrary constant parameters <Xj(0) = (vj(0), Dj(0), 7j(0) 
We assume that the nonlinearity (3 satisfies for all integers i > 

(2.7) l/? W ( s )l < s^"^+ forO<s<l 

(2.8) 



where p > 2 



<7 < 1 



|/?W(s)| < fors>l 
-. The main results of our paper are the following theorems. 



Theorem 2.1. Let (3 be as in (1.19), (1-20) with a sufficiently small value of parameter > 0. In 
addition, assume that (3{s 2 ) is a smooth function of s and that p + r > 1 + ^. Let 

N 

i>o( x ) = ^2wj(0,x) +R (x), 

3=1 

see (2.2) be initial data for NLS that satisfy the separation condition (2.17). Moreover, assume that Wj 
are defined in terms of ground states (positive radial variational solutions) of ENLS. Finally, suppose 
that for all a € R Ar ( 2n+2 ) with \a — cr(0)\ < c, the linearized operators Hj(cr) from (2.24) have zero as 
their only eigenvalue, no resonances, and no embedded eigenvalues in their continuous spectrum. Then 
there exists a positive e such that for Rq satisfying the smallness assumption 



(2.9) 



k=0 



for some integer s > S, there exists an admissible path a(t) with the limiting value a°° such that 



N 



4>(t,x) 



^2wj{t,x;aj{t)) 

3=1 



< 



(i+ty 



as t — > oo. Moreover, there exists uo £ L 2 so that 

N 



^2wj{t,x;a c 

3=1 



L 2 



as t 



oo. 



Remark 2.2. We assume that The smallness of the parameter 6 in the nonlinearity (3 depends 

on the values of <r(0). 

Our methods also allow treating initial conditions that are defined in terms of nonlinear eigenfunc- 
tions which are not ground states as well as more general nonlinearities (3. However, in contrast to 
the case of a ground state and the nonlinearities (1.19), in this case we cannot verify the convexity 
conditions as well as various conditions related to the spectrum of the linearized operators Hj{a). 
Therefore, we need to include them into the hypotheses. 
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Theorem 2.3. Let (3 be as in (2.7) and (2.8). Impose the separation and convexity condition see (2.17) 
and (2.16), the spectral assumption from Definition 2.4- Suppose ip is the solution of (2.1) with initial 
condition (2.2) where the Wj are generated by nonlinear eigenf unctions of ENLS (without assuming 
that they are ground states). Then there exists a positive e such that for Rq satisfying the smallness 
assumption (2.9) the conclusion of the previous theorem holds. 

In the remainder of this section we shall discuss the theorem and its assumptions in more detail. 
The soliton profiles Wj(t,x;a(t)): In addition to the one-soliton solutions Wj(t,x;o-j) with 
constant Oj introduced above, we need functions 

(2.10) Wj(t,x;aj(t)) = e ie ^' x ' u ^ (j){x - Xj (t; a(t)), aj (t)). 

The phase 9j(t, x; a(t)) and the path xj(t; a(t)) are defined in terms of the time-dependent parameters 
o-j(t) = (vj(t),Dj(t),jj(t),aj(t)) as follows: 

(2.11) e,(t,x;a(t)) = Vj (t)-x- J^\(\v \ 2 - a^s) ds + l3 (t) 

(2.12) Xj(t;a(t)) = f Vj (s) ds + D,(t). 

Jo 

Henceforth, the functions Wj(t,x;a(t)) or simply Wj(a(t)), correspond to the soliton moving along the 
time-dependent curve a(t) in the parameter space according to (2.11) and (2.12), while Wj(a) is the 
true soliton moving along the straight line determined by an arbitrary constant a as in (2.4) and (2.5). 

Admissible paths cr(t): We collect the individual parameter curves o~j(t) from above into a 
single curve a(t) := (<n(t), . . -,cr N (t)) C M.( 2n + 2 ) N . Given the initial value a(0) we introduce the set 
of admissible curves a(t) as those C 1 curves that remain in a small neighborhood of cr(0) for all times 
and converge to their final value cr(oo) = lim^ +oc a(t). We shall also impose the condition that for 
an admissible curve a(t) 

roo roc poo poo 

(2.13) / / \v j(t) ■ vj(t) — a(T)ja>j(r)\dT ds < oc, / / \vj (t)\ dr ds < oo 

JO Js JO Js 

for all 1 < j < N. Given an admissible curve a(t) we define the constant vector a°° in the following 
fashion: 

poo poo 

(2.14) v]°=Vj(oo), D°° = Dj(oo)- / vj^drds, 

Jo Js 

poo poo 

(2.15) 7 °° = 7i (oo) + / / ( Vj (T) ■ Vj (T) - aj(T)aj(T)) dr ds, af = aj (oo). 

JO Js 

Convexity condition: We impose the convexity condition 

(2.16) (d a c/)(-,a),^(-,a)) > 0, Va : min |a-a,-(0)|<c 

for some positive constant c. As noted in the Introduction, the convexity condition is closely connected 
with the issue of orbital stability of the individual solitons Wj(t, x; <7j(0)). 
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Separation conditions: Our theorem handles the case of so-called weakly interacting solitons. 
This means that the initial positions Dj(0) and initial velocities Vj(0) are such that for all t > one 
has the (physical) separation condition 

(2.17) |Dj(0) + Vj(0)t - D t (0) -v e {0)t\ > L + ct, Vj + £ = 1, . . . , N 

with some sufficiently large constant L and a positive constant c. Another assumption under which 
our theorems hold equally well is the condition of large relative velocities of the solitons. This means 
that 

(2.18) min \ vj — V(\ > L 
for some large L. 

Let a m i n = mini<j<fc ctj(0) — c. It will be understood henceforth that 

(2.19) a min L > | log e|. 

The small constant e appears in our theorem as a measure of smallness of the initial perturbation Rq. 

Spectral assumptions: We will write Wj(a(t)) = Wj(t, x;a(t)) = e^WH </>j(a(t)) where 
4>j(a(t)) = (j>(t,x;<Tj(i)) = (f)(x — Xj(t);aj(t)). Linearizing the equation (2.1) around the state w = 
YljLi w i-> ^ = w + R one obtains the following system of equations for Z = (^): 

(2.20) id t Z + H(t,a{t))Z = F. 
Here H(t,a(t)) is the time-dependent matrix Hamiltonian 

(2.21) H(t,a(t))=H + 

P(\wMm 2 ) + P\\wMt))\ 2 )\wMt))\ 2 a'(\wMt))\ 2 WM(t)) 



N 



P\\wMt))\ 2 )w 2 3 {a{t)) -f3{\ Wj {a{t))\ 2 ) ~ P'(\wMt))\ 2 )\™MtW 

N 



(2.22) =H + Y / V j {t,x-Mt)l 

with complex matrix time-dependent potentials Vj(t,x;a(t)) dependent on Wj and cr(t). The right- 
hand side F in (2.20) depends on &, w, and nonlinearly on Z. For a given constant parameter vector a 
we shall introduce the Hamiltonian H(t,a) 



(2.23) H{t,a) = H + 

0(\ Wj (o-)f) +(3>(\w j (a)f)\w j (o-)\ 2 (3>(\ Wj (a)f)w](o-) 

-(5'{\ Wj {a)?)w](a) -^{a)?) - 0'(\ Wj (a)\ 2 )\ Wj (a)\ 2 



E 



We refer to Hamiltonians of the form (2.23) as matrix charge transfer Hamiltonians. They are discussed 
in more detail in Section 12, as well as in [RSS]. Recall that Wj(a) denotes the soliton moving along the 
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straight line determined by the constant parameters Uj . The proof of our theorem relies on dispersive 
estimates for matrix charge transfer Hamiltonians that were obtained in [RSS], see also Section 12 
below. For these estimates to hold, one needs to impose certain spectral conditions on the stationary 
Hamiltonians 
(2.24) 

' \A - £ +/3(<Ma) 2 ) +(3>(Ma) 2 )<p J (o-) 2 ptf \ 



Hj(a) :-- 



where 4>j(o~) = (ft(x,aj), see (2.6). These Hamiltonians arise from the matrix charge transfer problem 
by applying a Galilei transform to the j th matrix potential in (2.23) so that this potential becomes 

2 

stationary (strictly speaking, this also requires a modulation which leads to the spectral shift ^- 
in (2.24)). We impose the spectral assumption as described by the following definition. 

Definition 2.4. We say that the spectral assumption holds, provided for all a € M. N ( 2n + 2 ) with \a — 
o"(0)| < c one has 

• is the only point of the discrete spectrum of Hj(a) and the dimension of the corresponding root 
space is 2n + 2, 

• each of the Hj(o) is admissible in the sense of Definition 12.1 below and the stability condition 

sup||e^ (<T) P s || 2 -,2 < oo, 

t 

see (12.3), holds (here P s is the projection onto the scattering states associated with Hj, see (12.2) ). 

While the second condition is known to hold generically in an appropriate sense, see Section 12, the 
first condition is more restrictive and not believed to hold generically. 

3 Reduction to the matrix charge transfer model 

For the sake of simplicity we consider the case of two solitons, i.e., N = 2. Setting 

(3.1) wi(<r(t)) +w 2 {a(t)) =w, tp = w + R, 
where Wj are as in (2.10), (2.11), and (2.12), we derive from (2.1) that 

(3.2) id t R + ^Ai? + (l3{\w\ 2 ) + (3'{\w\ 2 )\w\ 2 )R + (3'{\w\ 2 )w 2 R 

= -(id t w + ^Aw + I3{\w\ 2 )w) + 0(\w\ p - 2 \R\ 2 ) + 0(\R\ P ), 

using the assumptions (2.7), (2.8) on the nonlinearity (3. Observe that 

-, 2 

id t w + - Aw + I3{\w\ 2 )w = Yl [&■(*) ' x + 7j(*)W*(*)) + ie i0 ^ a{t)) V^(a(t)) • D 3 {t) 

3=1 

(3.3) -ie^( ff W)a a ^-(a(t))dj(t)l +0( Wl w 2 ). 
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In view of (3.1) one has 

id t R + ^AR + ((3(\w\ 2 ) + p'{\w\ 2 )\w\ 2 )R + p\\w\ 2 )w 2 R 

2 2 

= id t i? + iA^ + ^[/3(K.(a(^ 

3=1 3=1 

+0{ Wl (a(t))w 2 (a(t)))R. 
Rewriting the equation (3.2) as a system for Z = (R, R) therefore leads to 

(3.4) id t Z + H{a{t))Z = tW{a{t)) + 0{w lW2 )Z + 0{wiw 2 ) + 0{\w\ p - 2 \Z\ 2 ) + 0{\Z\ p ). 
Here H(a(t)) is the time-dependent matrix Hamiltonian from (2.21) and 

(3.5) tW(a(t)) = 
where 



/_ 



(3.6) / = £ [(vj(t) ■ x + j,(t))wMt)) + ie ie ^V4>Mt)) ■ Dj(t) - d a </>Mt))*i(tj 

3=1 

We shall assume that a(t) is an admissible path with initial values <r(0) in the sense of (2.13). Given 
an admissible curve a(t) we introduce the reference Hamiltonian H(t,a°°) "at infinity" 

H(t,a°°) = H + 

,0* ^^(K-(a-)| 2 ) + /3'(K-(0| 2 )K(0| 2 f3>(\ Wj (*™)\ 2 )w 2 (a™) \ 

[ } -P'(\ Wj (a^)\ 2 )w 2 (a^) -P{\w 3 {o-°°)\ 2 ) - /?'(|u>j(OI 2 )K-(OI 2 ) 

where a°° = (erf 3 , . . . , cr^ ) , a°° = (v°° , D°° , 7°° , a°° ) is the constant vector determined by the curve 
a(t) as in (2.14) and (2.15): 



v? = v 



poo poo 

(00), D°° = Dj(oo)- / Vj(T)dTds, 

Jo Js 

poo poo 

7 9° = 7i (oo) + / / (vj(T) ■ Vj(t) - dj(r)aj(T)) air ds, af = a^oo). 

JO Js 

Recall that Wj(a°°) is the soliton moving along the straight line determined by the constant parameters 

a™ . For j = 1, . . . , N we introduce the Hamiltonians 

(3.8) 

it (f n oo, _ R J /3(k J (^)| 2 ) + /?'(|^(a-)| 2 )|^(a-)| 2 (3>(\ Wj (o-°°)\ 2 )w 2 (c-°°) 

n j( i,a ) n +y ^^.^^^ -p{\ W j{o"°)\ 2 ) - /?'(K(OI 2 )KK°)| 2 

together with their stationary counterparts Hj(a°°) as in (2.24) with a = a°°. 

The following lemma relates the evolutions corresponding to the Hamiltonians Hj(t,a) and Hj(a) 
for an arbitrary a by means of a modified Galilean transformation. 
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Lemma 3.1. Let Uj(t,a) be the solution operator of the equation 

(3.9) id t Uj{t, a) + Hj(t, a)Uj(t, a) = 0, 

U j (0,a)=I 

and e ltH i( a } be the corresponding propagator for the time-independent matrix Hamiltonian Hj(o). Then 

(3.10) U 3 (t,a) = e; D] (t)^^)e ,ifl »^-(0,^,, Cj (0) 
where Gvj,Dj{t) is the diagonal matrix Galilean transformation 

and 

(3.12) Mj(t,a) = 

Proof. By definition, 



h) \Q Vj ,DAt)hJ 



e -l^-t-l{Vj-Dj+ 1: j) Q 

„,2 



iUj = iQ* Vj , Dj {t)M*{ty tH ^) M ^ 

(3.13) -a;,n.(t)A!K*)^(^)^ (<T) ^(0)^, ) D,(0). 



Clearly, 

iM*(t) 
whereas one checks that 



— rf e- 2 



i^-t+iivj-Dj+jj) q 



a i e-i-j-t-Uvj-Dj+yj) 



^ D] h) \ (v]/2 + v rP ) Q ^ D J 2/ 



Finally, we need to move Hj(cr) to the left in (3.13). We consider the differential operator separately 
from the matrix potential, i.e., 

/ Uj(x - xj(t)) e™*rtWj(x - xj(t)) 

j{h > ~ 0+ \-e- 2iB ^Wj{x-Xj{t)) -Ujix-x.it)) 

(3.14) H M ) = H 0+ (-J + % 

where Xj it),9j(t) are as in (2.5), (2.4), and Uj = P(4>j{a) 2 ) + /3 / (<^(<j) 2 )<^(c7) 2 , Wj = /3'(<^(<7) 2 )<£ 2 (cr). 
Note that, on the one hand, Mj commutes with all matrices in (3.14) that do not involve Uj, Wj. On 
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the other hand, one has 



H + 



2 < 



(t)fl 



i 

+ 2 



Z A 

e^'A^ 1 '"' e-ivjiDj+tvj) _ tv . _ _ j\t e ix- Vj e -i Vj iD J+ tv,) A/i(x - t Vj - Dj) 



-e~ l 2 * A 



^ e -«-^ f; ir, ■(/;,-/,;,) j 2 ( x _ ^ _ D .^J + e -i\t e -ix- Vj jvyiPj+tOj) A/ 2 (x - tVj - Dj)j 



1 („,2 i ^,2\ „* 



+ «j) <-,D,(*)/l " Vj-p9* Vi ,DMfl 



y-i(vj + aj) Q* v ., D .{t)f2-v r pQ* Vj:D .{t)h) 

Finally, we need to deal with the matrix potentials. Write M(t) := Mj(t,a) 
and set p = t\i>j\ 2 + 2x • Vj. Then (omitting the index j for simplicity) 



-iui(t)/2 







iw{t)/2 



M{t)Q^ D {t) 



U{- -vt-D) e 2ty W{- -vt-D) \ (fi 
-e~ 2ie W{- -vt-D) -[/(• -vt-D) ) \f 2 



e -iu>(t)/2 o 







HdV)U(- -tit- D)h + QvMfy ' w (- -M- D)f: 



Sv,D{t)e™W{- -tit- D)h - Q^, D {t)U{- -tit- D)j 2 



U^ D {t){e- iu >^/ 2 h) + We~ i( - V 2t +^e i ^ t ' + ^ +D ^\^ D {t)e i ^ t )/ 2 f 2 



_ We i(vH+2x-v) e i(u-20(t,-+tv+D)) 0t?jD (t)(e-^W/ 2 /i) - U ^, D {t)e^ t )/ 2 f 2 
U e i(20(t,-+tv+D)-w-p) W 
_ e -i(20(t,-+tv+D)-w~p) W _jj 



e -iw(t)/2 q 
e M*)/ 2 



9v{t)fi 

)v,D(t)fc 



Now 26 (t, - + tti+D)-p-LO = 2v-x + {\v\ 2 + a 2 )t + 2^ + 2v-D-t\v\ 2 -2x -v - u = by definition 
of to, i.e., uj = a 2 t + 27 + 2v • D. Adding these expressions shows that 

iUj(t, a) + Hj(t, a)Uj(t, a) = 0, 



as claimed. 



□ 



4 The root spaces of Hj(a) and H*(a) 

In view of Section 12 below (see in particular Definition 12.1 as well as (12.2)) we will need to understand 
the generalized eigenspaces of the stationary operators Hj(a) from (14.3). By our spectral assumption, 
see Definition 2.4 above, only generalized eigenspaces at are allowed. We denote these spaces by Nj (a) 
and refer to them as root spaces. Thus, Mj{a) = ker (Hj(a) 2 ) and by (12.2) one has the direct (but 
not orthogonal) decomposition 

L 2 (M 3 ) x L 2 {M. 3 )=M*{a) ± +M j (a), 

where M*(a) = kei (H* (a) 2 ) . The (nonorthogonal) projection onto M*(a) ± associated with this de- 
composition is denoted by Pj(cr). While the evolution e ltH ^ a ^ is unbounded on L 2 as t — > 00, it is 
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known in many cases that it remains bounded on Ran(P,-(<7)). In Section 12 this is referred to as the 
linear stability assumption. 

Proposition 4.1. Impose the hypotheses of Theorem 2.1 and let Hj(a) be as in (2.24). Then 

• The nullspace M* (a) ofHUa) is given by the following vector valued 2n+2 functions £J™ (x; a) , m = 



l,...,2n + 2: 



uf(x;a) 



u){x;a) = 4>j{x;a), 

2 2 
uAx;a) = i — d a (j)j(x;a), 

uj{x;a) = id Xm _ 2 4>j(x;a), 
uf(x;a) = x m ^ n - 2 4>j(x;a) 



H*(a)^(-;a) = 0, m = 3,.,n + 2, 

H*(a)^(.;a) = -2i^' n (.;a), m = n + 3, .., 2n + 2 



Lei 



J 



1 
-1 



Then J is an isomorphism between the nullspaces ofH*(a) and Hj(a). In particular, the nullspace 
of Hj(a) has a basis {J£™(-;cr) | 1 < m < 2n + 2}. Moreover, M*{a) is spanned by 



Jd ar Wj(t,x;a) for 1 < r < 2n + 2, 



where WAt,x-a) = ( w jf; x ' a \). 
• One has the linear stability property 



sup 

t 



< oo 



2^2 



where Pj(<r) is the projection onto M*(a) L as introduced above. 



Proof. 



□ 



For the case of monomial, subcritical nonlinearities these results go back to Weinstein's work on 
modulational stability [Wei]. 



5 Estimates for the linearized problem 

In (3.4) we obtained the system 
(5.1) 

id t Z+H(t,o-°°)Z = (F(CT(t))-i7(t,a^))Z+SI^(a(t))+0(wiw 2 )^+0(wiW2)+0(k| p ~ 2 |Z| 2 )+0(|Z| p ), 

The point of rewriting (3.4) in this form is to be able to use the dispersive estimates that were obtained 
in [RSS] for (perturbed) matrix charge transfer Hamiltonians, see also Sections 12 and Section 8 
in [RSS]. 
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Theorem 5.1. Let Z(t,x) solve the equation 

(5.2) id t Z + H(t,a)Z = F, 

Z(0,-) = Z (-) 

where the matrix charge transfer Hamiltonian H(t,a) satisfies the conditions of Definition 12.2. As- 
sume that Z satisfies 

(5.3) \\(ld-P 3 {a))M 3 (a,t)g v ^ D] {t)Z{t,-)\\ L 2 < 5(1 +<)"?, Vj = l,...,fc, 

u>rf/i some positive constant B, where Mj(a,t) and G Vj ,Dj(t) are as in Lemma 3.1. Then Z verifies the 
following decay estimate 

(5-4) \\Z{t)\\v + L°° < (1 + t)^ (||Z || L i nLa + \\F\\ + B) 

for t > with 

~ t \\F(s)\\ L1 ds + (l + t)% +1 \\F(t)\\ L * . 
In addition, we also have the L? estimate 

(5-5) l|Z(i)|| L2 < HZolUinia + lFl + 5 

For the proof see [RSS] and Section 12 below. In particular, note that (5.3) is related to the 
characterization of scattering states in Definition 12.4. 

In the applications the inhomogeneous term F is a nonlinear expression which depends on Z. 
Therefore, in addition to the estimates (5.4) and (5.5) we shall need corresponding estimates for the 
derivatives of Z. 

For an integer s > we define Banach spaces X s and y s of functions of (t, x) 



\F\ := sup 

t>0 k/o 



s 

(5.6) = sup(|hKt,-)||ff. + (l + t)S ^||VV(V)||l 2+ loo) 

(5.7) ||F||y. = supj]( r||V^(r,-)llL 1 dr + (l + t)f +1 ||V fc F(t,.)|| i2 ) 



t>o" Wo 

The generalization of the estimates of Theorem 5.1 is given by the following theorem (see Section 9, 
in particular Proposition 9.3 in [RSS] for the proof). 

Theorem 5.2. Under assumptions of Theorem 5.1 we have that for any integer s > 

s 

(5-8) \\Z\\ X . < II V*Z(0, OlUin^ + \\F\\y, + B 

k=0 
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We apply Theorem 5.2 to the equation (5.1). This will, in particular, lead to our main result, i.e., 
that HZ^Iloo < t~2 as t — > oo. We need to ensure that Z is a scattering solution relative to each of 
the channels of the charge transfer Hamiltonian H(t, a), in the sense of the estimate (5.3). Analogous 
to Buslaev, Perelman [BP1] this will be accomplished by an appropriate choice of the path cr(t), to be 
made in the following section. 

In order to prove existence of solutions Z and a we require another version of Theorem 5.2, which 
follows easily from that result. 

Remark 5.3. It is easy to see that a time-localized version of the previous theorem also holds. Indeed, 
let 



s 

(5-9) U\\ Xs (T) = SUp (U{t r )\\ H s +{l+tf^\\^ k ^-)\\L^L^ 

0<t<T V ^ > 

(5.10) \\F\\y s{T) = sup W /V fe ^,-)llL 1 rfr + (l + t)f+ 1 ||V fc F(t,.)|| L2 ). 

Then assuming (5.3) for < t < T with a constant Bt, one has 

s 

(5.H) \\z\\ MT ) < £ ll v * z ( . -)hw + \\F\\y,<T) + b t 

k=0 

Corollary 5.4. Let £™(i, x), 1 < m < 2n + 2, 1 < j < N, be a collection of smooth functions such that 

(5.12) sn P \\M 3 (a,t)g Vj , Dj (t)^ l (t,x)-^(-,a)\\ LlnL 2 < 5 

for some small 5 > and some given a. Let Z be a solution of 

(5.13) id t Z + H(t,a)Z = V(t,x)Z + F, 

(5.14) (Z(t),£™(t,-))=0 

for all t > 0, where V(t,x) is a smooth function that satisfies sup| 7 |<<, ||$?V(£, Ollz^nL 00 < S(l + i) -1 , 
with a nonnegative integer s for all t > 0. Then 

s 

(5-15) ||Z|U s <^||V fc Z(0,.)|| 
Proof. By (5.14), 

\{Mi{a,t)g Vj , Dj {t)Z{t),^{--,a))\ < \(M J (a,t)g v ^ Dj (t)Z(t),-M J (a,t)g V] ^(t)ij l (t,-) + 

<S\\Z(t)\\ 2+00 . 

Therefore, for all < t < T, 

\{Mj(a,t)g VjtDj (t)Z(t),^(t,--,a))\<S(l + t)-^ sup (1 + r)f \\Z(r)\\ 2+00 

0<t<T 

< 5(1 + i) _t \\Z\\x,<T) =■ (1 + Bt- 
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Hence, using (5.11) one sees that 

s 

\\ z \\x s (T) ^ ^ H^ 2 ^ ' Oll^nL 2 + \\ F \\y s (T) + ll^lly s (r) + b t 

k=0 

s 

<^2\\V k Z(0 r )\\ LlnL2 + \\F\\y s(T) + 5\\Z\\ X S (T), 
and the desired conclusion follows. □ 

6 Modulation equations 

In their analysis of the stability relative to one soliton, Buslaev and Perelman [BP1], [BP2], and 
Cuccagna [Cu] derive the equations for a by imposing an orthogonality condition on the perturbation Z 
for all times. More precisely, they make the ansatz 

(6.1) V = e ie ^ a{t)) {w{a{t)) + R) 

where e l9 ^' a<yt ^ w(a(t)) is a single soliton evolving along a nonlinear set of parameters. The removal of 
the phase from the perturbation R leads to an equation which is simply the translation of the equation 
involving the stationary Hamiltonian (2.24) to the point vt + D. This in turn makes it very easy to 
formulate the orthogonality conditions: At time t, the function R(- + vt + D) in (6.1) needs to be 
perpendicular to all elements of the generalized eigenspaces of all Hj(a)* as in (2.24), where a is equal 
to the parameters a{t) at time t. 

In the multi-soliton case the removal of the phases by means of this ansatz is not available, since 
distinct solitons carry distinct phases. As already indicated above, we work with the representation 

N 

ij(t)=Y,w j (t,a(t))+R, 

3=1 

which forces us to formulate the orthogonality condition in terms of a set of functions that is moving 
along with the Wj(t, a(t)). We now define these functions. 

Definition 6.1. Let a{t) be an admissible path and define 9j(t,x;a(t)) and Xj(t;a(t)) as in (2.11) 
and (2.12). Also, set (pj(t,x;a(t)) = cj)(x — xj(t; a(t)); aj(t)). Then we let 

v { ' ' ( } > \uf(t,x;a(t)) 



with 



(6.2) u)(t,x;a{t)) = Wj (t,x;a(t)) = e^ (i ' x;,j(i)) ^-(t, x; fj(t)) 
u 2 Jt, x- ait)) = — e ^(*.*;*(*)) d a 4>j(t, x- a(t)) 

Ctj 

uf(t,x;a{t)) = i e Wj(tw(t)) Xm _ 2 0j(t,z;<7(t)) for 3 < m < n + 2 

uf{t,x;a(t)) =e^ (i ' x;CT( * )) (x m - n - 2 -x™- n - 2 (t;a(t)))0 i (t,x;a(t)), for n + 3 < m < 2n + 2. 
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The following proposition should be thought of as a time-dependent version of Proposition 4.1. 
More precisely, if a is a fixed set of parameters, then one can define an alternate set of vectors, say, 
by applying appropriate Galilean transforms to the stationary vectors in Proposition 4.1. For example, 
take some £™ so that H^(g)£,™ = 0. Then the corresponding £™ satisfies 

with Hj(t,a) as in (3.8). Naturally, one would therefore expect that 

idtC? + H(a(t))iT = 0(&j) + 0(e~ ct ), 

where H(a(t)) is as in (2.21) (the exponentially decaying term appears because of interactions between 
solitons). The following proposition shows that this indeed holds, but as in [Cu] we will work with a 
modified set of parameters aj(t) = (vj(t),Dj(t),aj(t),jj(t)) where 



n 

(6.3) fat) = fat) + vf{t)x™{t,a(t)). 

771=1 

The point of this modification is that the SW(a(i)) term in (5.1) and (3.4) can be rewritten as 

k 

tW(a(t)) = Y,[fat)J^(t,x-,a(t))-^a :j (t)J^(t,x-a(t))] + 

k n 

(6.4) E [i>?(t)Jep 2 {t t xMt)) + -v^(t)J^ +n+ \t,x;a(t)) , 

j = l 777=1 

where £™ are as in Definition 6.1. This is of course due to the fact that passing to 7j allows us to 
change from x to x — Xj(t; cr(t)) in (3.6). 

Proposition 6.2. Let a(t) be an admissible path and define x; cr(t)) as in Definition 6.1. Then 

(6.5) id^+H^am] = 0(H\fa + \Dfa)) 

(6.6) ^ 2 + #;Hi))£f = it) + 0{d{\fa + \Dfa + \D 2 fa)) 

{Q.7)id t ty + H*{a{t))ty = 0(a(|^-| + |D^-| + |D 2 ^-|)) for 3<m<n + 2 

(6.8)^f + H*{a{t))£™ = -2i£™- n + 0@{\<f>j\ + \Dfa + |^ 2 ^|)) for n + 3 < m < 2n + 2. 

Here D refers to either spatial derivatives d x i or derivatives d a . Moreover, as in Definition 6.1, the 
function <pj needs to be evaluated at x — Xj(t;a(t)), ctj(t). 

Proof. This is verified by direct differentiation of the functions in Definition 6.1. □ 

The following proposition collects the modulation equations for the path a(t) that are obtained by 
taking scalar products of (2.21) with the functions £™ from Definition 6.1. This will of course use (6.4). 
The modulation equations are derived from the orthogonality assumptions, see (6.9) below. Observe 
that these assumptions need not be satisfied at t = 0. Nevertheless, as in Buslaev and Perelman [BP1], 
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one shows by means of the implicit function theorem that one can replace the initial decomposition (2.2) 
by a nearby one which does satisfy the orthogonality condition. This uses the smallness of the initial 
perturbation Rq, as well as the separation conditions (2.17) or (2.18). The details can be found in 
Section 11, see Lemma 11.1. In that section it is also shown that, conversely, given the modulation of 
Proposition 6.3 the orthogonality condition will propagate if satisfied initially. 

Proposition 6.3. Let Z satisfy the system (3.4). Suppose that for all t > 0, 

(6.9) (Z(t),^(t, •; a(t))) = for all j, m 

where £™ is as in Definition 6.1. Then the path a(t) := (vj(t), Dj(t),^j(t),aj(t)), j = l,..,n satisfies 
the following system of equations with matrix potentials V r (t,x;a(t)) as in (2.22): 

-2^(t)(^(a(t)), 9 Q ^(a(t)))+0(&||Z(t)|| L2+L co) = Xj(^ r (t,-;a(t))Z, $(t,-,a(t))) + 

(6.10) {[o{ Wl w 2 )z + o( Wl w 2 ) + o(\ w r 2 \z\ 2 ) + o{\z\ p )), •; 

2^(t)(^(a(t)), Q ^(t)))+O(*||Z(t)||^ +L =o) = J2(v r (t,-;a(t))Z, $(t, ■;*$))) + 

([0{w lW2 )Z + 0( Wl w 2 ) + 0(\w\p~ 2 \Z\ 2 ) + o{\z\ p )), $(t, •; <r(t)j), 

((0{w lW2 )Z + 0(w lW2 ) + 0(\w\ p ' 2 \Z\ 2 ) + 0{\Z\P)) , if +2 (t, ■; a(t)) 
Z?-(t)||^(,r(t))||l+0(5 r ||^(t)|U 2+J ,oo) = ^<T^(t, ^ + ^ +2 (t, -;^(t))> + 

( (o( WlW2 )z + o{ WlW2 ) + o(t«;r 2 izi 2 ) + o(i^n) , e; +m+2 (^, •; *(*))) 

Proof. Differentiating (6.9) yields 

<i$Z,^(t, ■;*(*))> = (Z,id t ^(t,-;a(t))). 
Taking scalar products of (3.4) thus leads to 

{Z,id t tiT)+{Z,H*{a{mT) = (tW{a{tV 

In view of the explicit expressions (6.2) one has 

(jg(t, ■;*(*)), #(<,-;*(*))> = -2^-(<7(i)), d a <f>j((T(t))) 

(J^(t,-;a(t)),^(t,-,a(t))) = for m + 2 

(J^(t,-;a(t)),^(t,-,a(t))) = for m + 1 

(Jef +2 (t,-;a(t)), ef +n+2 (t,-;cx(t))) = -2*||^Ht))||| for 3 < m < n + 2. 

Therefore, the proposition follows by taking inner products in (6.4). Note that the terms containing 
cr||Z(t)|| L 2 +L oo appear from Proposition 6.2. □ 
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7 Bootstrap assumptions 



The proof of our main theorem relies on the bootstrap assumptions on the admissible path a(t) and 
the size of the perturbation Z(t,x) = (^'^)- i n the norms of the spaces X s defined in (5.6). 
Bootstrap assumptions 

There exists a small constant 5 = 5(e) dependent on the size of the initial data Rq and the initial 
separation of the solitons Wj(0, x; <r(0)), see (2.17), and a sufficiently large constant Cq such that for 
some integer s > § 

(7.1) \a(t)\ < 5 2 (1 + t)~ n , Vt>0, 

(7-2) \\Z\\ X , < 5C, 1 

Remark 7.1. The bootstrap assumption (7.1) together with the definition (6.3) implies that 

(7.3) |7(t)| <S 2 (l + ty n+1 

Remark 7.2. The bootstrap assumption (7.2) together with Lemma 9.4 implies that 

(7.4) \\Z(t)\\L-<SCo\l + t)-S, 
(7-5) \\Z{t)\\ H s<5C^ 

The bootstrap assumption (7.1) strengthens the notion of the admissible path. In particular, it 
allows us to estimate the deviation between the path Xj(t;a(t)) corresponding to the path a(t) and 
the straight line Xj(t,a°°) determined by the constant parameter a°° which was defined from a(t) in 
(2.14) and (2.15). This estimate will play an important role in our analysis. 

Lemma 7.3. Let a(t) be an admissible path satisfying the bootstrap assumption (7.1) and let a 00 be a 
constant parameter vector as in (2.14) and (2.15). Then 

(7.6) \ Xj (t; a(t)) - tvf - Df\ < 5 2 {l + t)- n+2 
Proof. By our choice of v°° and D°° one has that 

roc roo roo 

\ Xj (t-o(t))-tvf -Df\ < / / \v j (r)\dT+ \Dj(s)\ds 

Jt Js Jt 

and the lemma follows from (7.1). □ 
We then have the following corollary. To formulate it, we need the localizing functions 

Xo(x) = exp ( - -a min (l + |x| 2 )5^ 

k 

(7.7) X(t,x;a°°) = ^ X o(x - xj(t; a 00 )). 

3=1 

Here a m \ n > satisfies inf t >o,i<j<fe oij(t) > a m \ n for any admissible path a(t) starting at do- The 
exponent a m - m arises because of the decay rate of the ground state of (2.6). 
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Corollary 7.4. Let a(t) be an admissible path satisfying the bootstrap assumption (7.1). With the 
parameters a°° as in (2.14) and (2.15) one has 



(7.8) H(t,<r°°)-H(*(t)) <5 2 (l + t) 2 - n X (t,x;o-n, 

where H(t,a°°) and H(a(t)) are the Hamiltonians from (3.7) and (2.21). 
Proof. The difference 

H(t,v°°)-H(v(t)) 

is a sum of matrix valued potentials that are exponentially localized around the solitons Wj(a(t)) 
or Wj(t, <7°°), respectively. By the previous lemma, we can assume that all the potentials are localized 
near the straight path x(t, a°°) = (xi(t; cr°°), . . . ,a;jv(t; Since t>°° = Vj(oo), a°° = <x,-(oo), 



(7.9) 



t 

[H(t,o-°°)-H(o-(t))}\<J2\tvr+D°°- [ v,(s)ds-D 3 (t) Xo (x - Xj (t; a 00 )) 

' .7 = 1 J ° 



(7.10) 



I ]_ /*00 ^ /*£ rCxD 

The term (7.9) arises as the difference of two paths, whereas (7.10) is the difference of the phases, i.e., 

uiOjfrx;* 00 ) _ e i0j(t,x;a(t))\ 

In view of the definitions of Dj,jj from (2.14) and (2.15) one has 
(7.11) 

h 

f*00 P OO /* OO 

H{t,o™)-H{a{t))\<jXj t I \v 3 {r)\dT + J^ \D 3 {s)\ ds) X o{x - Xj(t, a)) 

h 

, poo poo poo poo . 

+ J \vj{T)-v j {T)-a j {s)a j {s)\dTds + J \ij(s)\d8 + J \vj(s)\ ds\x\) X o(.x - Xj (t, a 00 )) 

<<5 2 (i + *) 2 -"x(M;0- 

For the final inequality one uses (7.3) and the fact that 

\x\xo(x - Xj(t;a°°)) < t. 

The corollary follows. □ 
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8 Solving the modulation equations 



Our goal is to show that the system in Proposition 6.3 has a solution a(t) that satisfies the bootstrap 
assumptions (7.1). This requires some care, as the right-hand side in Proposition 6.3 involves the 
perturbation Z. We will therefore first verify that the system of modulation equations is consistent 
with the bootstrap assumptions (7.1) and (7.2). In what follows, we will use both paths a(t) and a(t). 
By definition, see (6.3), 



ds. 



r°° r 1 , — - n 

W) = - / 7iW + ~ E vf(s)xf( S -a{s)) 

Jt m=l 

The integration is well-defined provided a satisfies the bootstrap assumption. Indeed, in that case 
^ (1 + *)~ n an d since \xj(t;a(t))\ <l + t, the integral is absolutely convergent. Finally, recall 
the property (7.3) of the derivatives. 

Lemma 8.1. Suppose the separation and convexity conditions hold, see (2.17) and (2.16). Leta,Z 
be any choice of functions that satisfy the bootstrap assumptions for sufficiently small 5 > 0. If the 
inhomogeneous terms of the system (6.10) are defined by means of these functions, then this system 
has a solution a that satisfies (7.1) with 5/2 for all times. 

Proof. By the nonlinear stability condition (2.16), the left-hand side of (6.10) is of the form Bj{t)aj{t) 
with an invertible matrix Bj(t). The O-term is a harmless perturbation of the matrix given by the 
main terms on the left-hand side, provided 5 is chosen sufficiently small. This easily follows from 
the smallness of Z given by (7.2). We need to verify that the right-hand side of (6.10) decays like 
<5 2 (1 + t)~ n . We consider only the first equation in (6.10), the others being the same. The terms 
(V r (t,a)Z,^j(t, -;a(t))) for r ^ j and w\w 2 are governed by the interaction of two different solitons. 
In view of the separation condition (2.17) and the exponential localization of the solitons, we have 

(8-1) |dy(t)| < e —n(L +ct)(1 + \\ Z (t)\\ L i +L „) + ||Z(t)|| 2 i2+Loo + \\Z(t)\f L2+Loo 

(8.2) KSC-^l+tr^e + SC' 1 + 6^0-^) < Q\l + t y n 

where we have used the estimate (7.4), the condition (2.19), La m \ n > \ loge|, and that p > 2. □ 

More generally, the estimates leading up to (8.1) also yield the following result. The proof is implicit 
in the preceding one and is therefore omitted. 

Lemma 8.2. Let a(t) be an admissible path satisfying the bootstrap assumption (7.1) and Z(t) be an 
arbitrary function in X s . Define the function X as a solution of the equation 

(8.3) {tW(a(t)), if{t, ■■o{t))) = (G(Z(t),a(t)),^(t, •; a(t))) + (Qf(t, •; a(t)), Z(t, •)}, 
where 

(8.4) G(Z(t),a(t)) = 0( Wl (a(t))w 2 (a(t))Z + 0{ Wl {a{t)w 2 {a{t)) + 0(\w(a(t))r 2 \Z\ 2 ) + 0(|Z| P ), 

(8.5) nf(t,x;a(t)) = 0(d(\4>j\ + \D^\ + \D 2 4>j\)) +Y,Vr(t,x;a(tmj l (t,x;a(t)). 
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Then 

(8.6) \t(t)\ < (1 + t)~"(^ 2 + C\\Z\\ 2 Xa + CIIZII^) 

Remark 8.3. The functions G(Z(t), cr(t)) and x; c(t)) arise as follows. In Section 11 we will 

rewrite the Z equation in the form 

id t Z + H(a(t))Z = t\V(a(t)) + G{Z(t),a(t)), 

The quantity £t™(t, x; cr(t)) is defined via the equation 

id t £™(t, ■; a(t)) + H*(a(t))^(t, ■; a(t)) = S?${t, •; a(t)) + ^(t, •; a(t)), 

where the matrix S collects the terms £™(i, •; cr(t)) on the right-hand sides of (6.5)-(6.8). However, the 
previous proof does not require the explicit form of G or QJ 1 . 

9 Solving the Z equation 

In this section we verify the bootstrap assumptions (7.2) for the perturbation Z. This together with the 
already verified bootstrap estimates for a will also lead to the existence of the function Z(t) asserted 
in our main result. At this point we recall the imposed orthogonality conditions (6.9) 

(9.1) (Z(t),^(t,-,a(t))) = for all j,m 

with x; <r(i)) is as in Definition 6.1. We next rewrite the equation (3.4) for Z in the form 

(9.2) idtZ + H(t,a°°)Z = F, 

(9.3) F = (H(t, O - H(a(t)))Z + tW(a(t)) + 0{ Wl w 2 )Z + 0{ Wl w 2 ) + 0{\w\ p ~ 2 \Z\ 2 ) + 0(\Z\ P ) 

with the reference Hamiltonian H(t, a°°) as defined in (3.7). To verify the bootstrap assumption (7.2) 
we need to apply the the dispersive estimate for the inhomogeneous charge transfer problem stated in 
Theorem 5.2. The following lemma shows that the orthogonality conditions (6.9) and the bootstrap 
assumptions (7.1) , (7.2) imply that Z is asymptotically orthogonal to the bound states of H*(a°°), as 
required in Theorem 5.2. 

Lemma 9.1. Let Z be an arbitrary function in X s satisfying the orthogonality conditions (9.1) with 
respect to an admissible path o~(t) obeying the bootstrap assumption (7.1), and so that Z verifies the 
bootstrap assumption (7.2). Then Z is asymptotically orthogonal to the null spaces of the Hamiltonians 
H*(a°°) in the sense of (5.3). In fact, 

(9.4) \\P Nj (o-°°)g vT , Dr (t)Z(t,-)\\ L i <«5 3 (l + t)-§-\ Vj = l,..,fc 

Proof. By the assumption Z(t) is orthogonal to the vectors £™(t, x; a(t)) introduced in Definition 6.1, 
while (9.4) is equivalent to the estimates 

\(9^, Dr (t)Zr(--,v°°),Z(t))\ <* 3 (l+t)-*-\ Vj,m 
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Here £™(x; cr°°), defined in Proposition 4.1, refer to the elements of the null spaces Mj{a ao ) of the sta- 
tionary Hamiltonians H*-{o°°). The desired estimate would then follow from the bootstrap assumption 
(7.2), in particular (7.4), and the inequality 

(9.5) \K T ,D T mr{-, o - er(t, -Mm* < ^ 2 a + r 1 

The vectors £™ are composed of the functions derived from the bound state <p. In particular, £j = (^) . 
Therefore, 

\g* vr!Dr (t)$(x;a°°) - $(t,x;a(t))\ =2\e^ v r^i(K\ 2 -\ a T\ 2 ) t ^r) 4>( x - v ft - Df)- 

iC- Vj (tyx~\}{\v 3 {T)\^- aj {rr)dr+ l3 (t)) 

(9.6) e o (f)(x - Xj (t,x;a(t)))\ 

According to Lemma 7.3, \xj(t, x; a(t)) - vft - Df\ < 5 2 (l + t)- n+2 . Similarly, (7.10) of Corollary 
7.4 gives the estimate for the difference of the phases appearing in (9.6) 

^ e iSj(t,x;a°°) _ (t,x;<r(t) | < _|_ 

The estimate (9.5) follows immediately since n > 3. □ 

We now in the position to apply Theorem 5.2 to establish the improved X s estimates for Z(t). 

Lemma 9.2. Let Z be a solution of the equation (9.2) satisfying the bootstrap assumption (7.2) with 
some sufficiently small constants 5 and Cq 1 . We also assume (due to Lemma 8.1) that the admissible 
path a(t) obeys the estimate (7.1). Then we have the following estimate 

(9-7) \\Z(t)\\ Xs < 5 -C Q l 

Proof. Perturbation Z is a solution of the inhomogeneous charge transfer problem (9.2) 
id t Z + H(t,a°°)Z = F, 

(9.8) F := (H(t, a°°) - H{a{t)))Z + t\V(a(t)) + 0( Wl w 2 )Z + 0( Wl w 2 ) + 0(\w\ p ~ 2 \Z\ 2 ) + 0(\Z\ P ) 

Lemma 9.1 shows that Z is asymptotically orthogonal (with the constant S 3 ) to the null spaces of the 
Hamiltonians H*{a°°). Therefore, Theorem 5.2 gives the estimate 

s 

(9-9) \\Z{-)\\x s < £ HV^oIIlW + \\F\\y. + S 3 

with 

(9.10) \\F\\y. = sup (V f ||V fc F(r, -)|| L1 dr + (1 + t)^\\F(t, .)\\ H .) 

By the assumptions on the initial data J2k=o \\^ k ZoWl^dl 2 < 6 << Therefore, to obtain the 
conclusion of Lemma 9.2 it would suffice to verify that 

(9.H) \\F\\y s <5 2 

with F defined as in (9.8). The estimate (9.11) relies on the following lemma and the bootstrap 
assumptions (7.2) on Z(t). □ 
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Lemma 9.3. Let a(t) be an arbitrary admissible path satisfying the bootstrap assumption (7.1) and 
Z{t) be an arbitrary function in X s . Then the nonlinear expression F defined in terms of the path a{t) 
and Z(t) as in (9.8) obeys the estimate 

(9-12) \\F\\y s < 5 2 + \\Z\\ 2 Xs + \\Z\f Xs 

9.1 Algebra estimates 

In this section we establish several simple lemmas designed to ease the task of estimating the y s norm 
of F = F(Z, w, a) in connection with the X s norms of Z. 

We start by formulating a version of the Sobolev estimate tailored to the use of the space L 2 + L°°. 

Lemma 9.4. Let s be a positive integer. Then for any nonnegative integer k < s and any q £ [2,%], 
where 

(9.13) - = \~—, *fk>s-^ 

q k 2 n 2 

q k = oo, if k <s - - 
and q € [2, oo) if k = s — ^ the following estimates hold true 

s 

(9.14) l|v fc /lk 9+ L- < £ \\v l f\\ L 2 +L ~ < (i + r ? ll/lk 

1=0 

In particular, if s > ^ 

(9.15) n/iiL~<(i+r%ik 



Proof. By duality and density it suffices to show that 

s—k 

\\f\\L^<Y,W vl fh^' 

1=0 

The L 1 estimate is trivial while the the estimate for the I? norm follows from the standard Sobolev 
embedding W k ~ l ' q ' C L 2 , which holds for the range of parameters (k, I, q) described in the Lemma. □ 

Next are the estimates of the nonlinear quantities arising in (9.8) in terms of the X s norm. 
Lemma 9.5. Let j(t) be a smooth function which obeys the estimates 

(9.16) |7 W (t)|<t<^-4+ 
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for some p > 2 + ^ and all non-negative integers £. Here r + = r if r > and r + = if r < 0. Then 
for any s > ^ and any non-negative integer k < s 

(9.17) l|V fc ( 7 (|/| 2 )/) H^i < (1 + t)-\\\f\\ p Xs + Wfft 1 ^- 1 ' 2 ^), 

(9.18) l|v fc ( 7 (|/! 2 )/) || L2 < (i + t)-%-H\\f\f Xa + ||/|lS +1+(p_1 " 2fe)+ ) 

In addition, if 7 is a smooth function obeying (9.16) /or some p > 2 and ("(x) is an exponentially 
localized smooth function, then for any q G [1,2] 

(9.19) nv fc (c7(i/i 2 )/)iiL, < (i+r n (ii/i^ + ii/n| +1+(p ^ 2fe)+ ). 

Remark 9.6. It will become clear from the proof below that if the function 7 satisfies (9.16) for some 
p > 2 + ^, then the estimate (9.17) holds with a better rate of decay in t. In particular, 



(9.20) ini vfe ( 7(|/|2)/ ) 

Proof. By Leibnitz's rule 



L dt< \\f\\ p Xs + 



V fe (7(l/| 2 )/) =E E ^7 W (l/| 2 )V mi /V m2 /...V m -+V 

<?=0 mi+...+rri2^-|-i=fc 

with some positive integer constants and non- negative vectors rh = (mi, . . . , m 2 ^ + i). We may 
assume that m^+i > m 2 e > ■ • • > mi. Define 

(9.21) a mr = (I - £_^)-i 

for m r > s — § and q mr = 00 otherwise. With the above definition the Sobolev embeddings 

(9.22) H s C W mr ' qmr , W s ' 2 + W s,co C W nLrSmr + W nLr,0 ° 
(recall that m r < k < s) 3 hold true by Lemma 9.4. Then 

I|v fc ( 7 (l/I 2 )/)||L 1 < ||7(!/l 2 )lv fe /r^ll LlnL7 - T l|vVll! 2+LOO + 

(9.23) J2 E ll7 (0 (l/| 2 )V m V... V m -/|| LinL9 ^ +i (||/|| i2+L oo + ||V/||l»+l-), 

where the final term inside the parentheses in (9.23) arises via the second embedding in (9.22). We 
claim for £ > that there exist two sets of parameters and g^ r for r = 1, . . . , 2£ such that 

(9.24) 2<q 1 J r <q mr , W = 1, . . . ,2£, 

21 21 
r=1 9m r r=1 9m r Qrn 2 e+i 



3 In the case of m r = s — § the value of q mr can be set arbitrarily large, but the embedding fails for q mr = 00. However, 
since the following argument has some "slack" , we can allow ourselves to still set q mr — 00 for simplicity. 
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To prove the claim we let r be the number of m r . for r = 1, . . . , 21 such that m r > s — S. Observe that 



2 • 

21 



E m r < k - m 2 e+i 



r=l: m r >s — ^ 



Therefore, 



2^ 

^q mr ~ 2 n ~ n 2 n 

" n 2> n 2-2-q' m2e+i 

The inequality in the second line above follows since m^+i > s — \ , which holds if r > 0. On the 
other hand, 

2£ 

it 2 " 

and the claim immediately follows, provided that t > 0. Thus using the sequence q\ Ur to handle the 
L 1 norm in (9.23) and for the L qm2l + l norm, we obtain 



|V fc (7(|/| 2 )/)lUi < ll7(l/| 2 )|VV| 1 ^ll LlnL7 - T l|VVll| +Loc + 



E E ii7 (£) (i/i 2 )iiL~ iiv-zii^^ . . . iiv-/n L ^ nA IIV/II^+lo 

£=1 mi+...H-m 2 ^+i=fc 



By Holder's inequality 



2^ 

f. 1 ' ' 11 ' 1 1 - HJ II ^ , -C I •/ ff." 



IWI/I 2 )^/! 1 --!!,^^ < llvVlt-ll/ir 1 - < im|P 



provided that p > 2 + ^, which is dictated by the condition that (p — 1)^^ > 2. Finally, using the 
property (9.24) together with the estimate (9.16) we obtain 

I|v*( 7 (I/I 2 )/)|Ili < ll/ll^ f liv fe /ll| +i oo + (ll/ll^ 1 + ll/ll^ +(p_1 " 2fc)+ )IIV7llL 2+ Loo 

^^(ii/ii^ + ii/iiS^" 1 " 2 ^)- 

Similarly, we estimate 

(9.25) ||V fc ( 7 (|/| 2 )/)|| L 2 < ||7(l/| 2 )llL2nL«>l|V fe /||L2 +L o 

k 



> + 



E E ii7 w (i/i 2 )v m v . . . v^/n (ii/Hl» + l- + ii v/Hl»+l-) 

<=lmi+...+m 2<+1 =lt L 2 nL*™27+i 2 
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To estimate the first term in (9.25) we note that 

ii7(i/i 2 )iu 2 nLoo < ii/iiej^d < ii/ii^ii/Hl^ < (i+r f(j, - 2) ii/ii^\ 

where the second inequality follows from interpolating L 2 ^ -1 ) between 1? and L°° and the last in- 
equality is a consequence of Lemma 9.4 and the definition of the space X s . Thus 

(9-26) ll7(l/| 2 )llL 2 n^l|V fe /llL 2+ L» < (1 + ||/||^. 

Furthermore, using definition (9.21) we have that 

Qm 2e+1 - 2 s- m 2 e+i ' 

Then 

(9.27) || 7 W(|/| 2 )V-V---V^/|| L2nLi= ^ < ||V m VllL- 1+ LocX 

||7 W (l/| 2 )V m2 /...V m -/ll 

L 2 nL a ~ m ' 2e +'t nL qm i~ 2 nL 9m i ( s - m 2i+i)- n 

Using the definition of q mi from (9.21) and the assumption that m\ < m 2 ^ +1 we infer that the last 

norm reduces to the one of the space 

9 — ^ — n 
L z n L a - m 2e+i , for mi < s - -, 

L 2 (~1 ^"-'«2f+i+('»—i-f) for mi > s 

2 

We now let r be the number of m r for r = 2, ..,2£ such that m r > s — Observe that since s > § 
and s > k 

21 

T ST — k + 771-2^+1 + mi 



r=2 



9m r 2 n 

s 1 fe — m 2m - mi 
= - r( n"2 ) + n 

^ . ,s-m 2 £+i (s - m 2 e+i - mi) + s - § 

< mm{ , -j 

n n 



On the other hand 



^2 2 - 2 

r=2 

It therefore follows that there exist 2 sets of parameters q^ and q^ for r = 1, . . . , 2£ such that 



(9.28) 2<q^ 2 r <q mr , Vr = 2,...,2£ 



21 

El (s - m 2 £ + i - mi) + s - § s-m 2 e+i 

-TT- = or . 

r=2 Tmr n n 
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In either case, with the help of Lemma 9.4, we can estimate 
||7 W (l/| 2 )V m2 /-..V m2 VII 



L 2 nL s ~'™ 2 *+ 1 nL 9m i~ 2 nL 9m i (s ~ mM + 1, ~ 



It therefore follows that the second term in (9.25) is 

k 



(9.29) 



<(i+r n E n/ii t 1_2 ' )++2 ' +1 - 



Now combining this with (9.26), and using the condition that p > 2 + -, we infer that 



|v fe ( 7 (l/l 2 )/)ll^ < + 



p x ■+\\f\\T 1+{p - 1 - 2k)+ )- 



The proof of (9.19) proceeds along the lines of the argument for the 1? estimate (9.18). We first 
observe that since C,ix) is an exponentially localized function, the L q estimate for 1 < q < 2 can be 
reduced to the L 2 estimate. We then note that the condition that p > 2 + ^ was only used in the 
estimate (9.26) which now takes the form 

||7(l/| 2 )CI|L 2 nLoo||V fc /l|L 2+ Loo < ||7(l.f| 2 )IM|V fc /|| L2+Loo 

<lllfl^ 1 || Loo ||V fe /llL 2+ L- 



□ 



The remaining estimates already have the desired form (9.29). 
9.2 L 1 estimates 

In this section and the following we prove Lemma 9.3. We start with the verification of 



s poo 

E / \W k F{t,-)\\ L ,dt<8 2 + \\Z\\% + \\Z\ 



with F as in (9.8). By Corollary 7.4 we have 



(9.30) 



H(t,v°°) - H(a(t)) < 5\l + tf- n X {t,x;o™), 



where x(t, x; <7°°) is a smooth cut-off function localized around the union of the paths Xj(t;a°°) = 
Vj°t + Dj°. Moreover, the spatial derivatives of the above difference also satisfy the same estimates. 
Using the bootstrap assumptions (7.2) we obtain 



(9.31) 



S /*00 . r ... S f'OO 

E / V\[H(T,a°°)-H{a{T))\Z{T)) dr < 6 2 E / l|V fc Z(r) \\ L * +L ~ (1 + r) 2 ^ dr 



<5 2 \\Z\ 



A',. 



/•oo 

I ( 



l + T) 2 - n -^dr <5 2 \\Z\ 



X s 
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The term T»W(a(t)) obeys the point wise bound 

\tW(a(t))\ <max|^)| X (t,x;0 

3 

This can be easily seen from the equation (6.4) and Lemma 7.3. The same estimate also holds for the 
spatial derivatives of the quantity above. Thus, with the help of the already verified estimate (7.1) we 
infer that 



(9.32) 



£ 

fc=0 



OO t* oo 

"v k {tW{a{r)))\ Li dr<8 2 (1 + r)~ n dr < 6 2 



The estimates for the 0(w\W2)Z and 0{w\W2) terms in (9.8) are straightforward due to the separation 
and the exponential localization of the solitons w\ and W2, e.g., 



(9.33) r\\0(v k (w lW2 ))\\ Ll < I 

k=0 Jo V / J 



Li ^ / e- a ^ L+c ^ dr < < < S 2 



Here we have used the separation assumption (2.17) and the condition (2.19), a m \ n L > | loge|. 

The exponential localization of the multi-soliton state w, the bootstrap assumptions (7.2) and the 
estimate (9.19) of Lemma 9.5 yield the estimate 



(9.34) 



S si oo /* oo 

\\oh k (\ W r 2 Z 2 ))\\ L ,dr<\\Z\\ 2 x l {l + T)- n dT<\\Z\\ 2 Xs . 
k=0 Jo v > Jo 



Finally, with the help of (9.17) (more specifically using the improvement (9.20) of Remark 9.6), we 
obtain 



(9.35) 



\\V k (Z p (T))\\ Li dr<\\Z\ 



9.3 L 2 estimates 

In this subsection we establish the estimate 

\\F(t, -)\\ H . < (1 + t)-?- 1 (5 2 + \\Zf Xs + \\Zf Xs ). 

The arguments follows closely those of the previous section. Using the estimates (9.30), (9.19) and the 
bootstrap assumptions (7.2) we obtain 



(h(t, O - H(a(r))) Z(f)\\ <S 2 (1 + t) 2 ~ n £ || V k Z(t, •) \\ L2+Loo 

k=0 

(9.36) < 5 2 \\Z\\ Xs (1 + t) 2 -"-f < 5 2 \\Z\\ Xs (1 + ty^ 1 
where the last inequality follows since n > 3. Similar to (9.32) 

(9.37) \\±W(a(t))\\H* < S 2 (l + t)~ n <5 2 (1 + t)'^ 1 
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The estimates for the 0{w\W2)Z and 0{w\w-i) terms again follow from the separation and the expo- 
nential localization of the solitons w\ and u>2, 

(9.38) \\0{wiw 2 )\\h' <e~ a ™^ L+c ^ <5 2 {l + t)-%-\ 

The exponential localization of the multi-soliton state w together with the estimate (9.19) of Lemma 
9.5 and the bootstrap assumption (7.2), also give the estimate 

(9-39) \\0(\wr 2 Z 2 )\\ H s < \\Zf Xs (l+t)- n < \\Zf Xs (l+t)-%-\ 

Finally, using the estimate (9.18) of Lemma 9.5, we obtain 

(9-40) \\ZV(t)\\ H s<(l + t)-^\\Z\\ Xs . 

This completes the proof of Lemma 9.3. 

10 Scattering 

In this section we intend to prove the last part of Theorem 2.1. More precisely, we shall show that 
there exists uq € L 2 so that 



N 



-> 

L 2 



as t — > oo. Observe that the equation (5.1) for Z = (R, R) and the Corollary 7.4 allow us to write the 
solution ip(t) in the form 

N 

ij)(t) = ^2w j {t,x-<j°°) + R{t,x)+Q{t,x), 

3=1 

where the function Q is spatially localized and decays in time with the rate t~ n+2 and thus ||Q(£)||i,2 — > 
as t — > oo. The function R verifies an inhomogeneous Schrodinger equation 

id t R + ^AR = V(t,x)R + F, R(0,x)=R (x) 

with spatially exponentially localized potential V. According to Lemma 9.2 the inhomogeneous term 
F satisfies the estimate 

(10-1) \\F{t,-)\\ L ,<{l+t)-^8 2 , 
while the solution R has been shown to obey the dispersive estimate 
(10.2) \\R(t,-)\\ L 2 +L ^<(l+t)-%5. 

The standard scattering theory argument show that the desired function no can be constructed as the 
limit 

Uo = R -i lim f e" i i A (V(s, -)R(s, •) + F(s, •)) ds. 
Jo 

The estimates (10.1) and (10.2) together with the localization of V guarantee the existence of the limit. 
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11 Existence 



In Lemmas 8.1 and 9.2 we established the estimates 

(11.1) |^|<^ 2 (l + i)-" 

for the admissible path a(t) and 

(H-2) \\Z\\x s <5 2 

for the solution Z(t,x) of the nonlinear inhomogeneous matrix charge transfer problem (9.2), under 
the bootstrap assumptions (7.1), (7.2) 



(ii.3) \z\ < 5 2 {i + ty 



(H-4) \\Z\\ X , < 5C, 1 



and the condition that Z is asymptotically orthogonal to the null spaces of the Hamiltonians Hj(a) 
with the constant J 3 . In this section we shall show that these are sufficient to establish the existence of 
the desired admissible path and the perturbation R. We prove existence by iteration. We shall define 
a sequence of admissible paths a^ n \t) and approximate solutions Z^ a \t) for n = 1, . . . according to 
the following rules. First, we write the Z equation (9.2) and (9.3) in the form 

(11.5) id t Z + H{a{t))Z = t\V{a{t)) + G{Z, a{t)) 

where G(Z,a(t)) = 0( Wl w 2 )Z + 0( Wl w 2 ) + 0(\w\p- 2 \Z\ 2 ) + 0(\Z\p). Set 

<j«(i) = <r(0), Z(°)=0 

where <r<n is to be understood as the constant path coinciding with the initial data <r(0) common to 
all admissible paths. We now define functions Z l (t, x) and cr 2 (t) to be a solution of the following linear 
system 

(11.6) 

idtZW + H(aW(t))ZW = £WW(<rW(t)) + G(Z«V (1) (t)), 
Z^\o,x) = Z (x), 
(11.7) 

£Ww(*i 1 \t)),&(t r] *( 1 \t))) = (G{zW 



Here f2™(i, x; a^\t)) is defined via the equation 
^(V^^+iHa^ 

where the matrix S collects the terms •; crW (t)) on the right-hand sides of (6.5)-(6.8). Thus, 

using Propositions 6.2 and 6.3, we have 

^^(OHO^^-K^ 
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Observe that (11.6) arises from the nonlinear equation (11.5) by replacing a(t) with the already defined 
path <r^(t) as well as Z on the right-hand side with = 0. The equation (11.7) determining a^ 2 \t) 
ensures that (zW {t),gj*(t, •; cr (1) (*))) = 0. Indeed, taking scalar products of (11.6) with £™(t, •; o-W (*)) 
and using Propositions 6.2 and 6.3 yields 

where E(i; a^\t)) is the vector of (£™(i, •; a^\t)), Z^\t, •)) and where S is the constant matrix defined 
above. As long as 

(11.8) S(0;a( 1 )(0))=S(0; < 7(0))=0 

one therefore has E(i; a^\t)) = for all t > 0. Generally speaking, (11.8) need not be satisfied. How- 
ever, using the fact that the initial perturbation Rq is small, we proceed as in [BP1] Proposition 1.3.1 
to show that one can modify the initial splitting in such a way that it does hold. More precisely, one 
has the following lemma. 

Lemma 11.1. Let ^(0) = YljLi w j(®i x i a '(ty) + -^o with Rq small as in (2.9) and assume either 
one of the separation conditions (2.17) orj^.18) as well as the convexity condition (2.16) in a small 
neighborhood U of c(0). Then there exist cr(0) € U such that in the decomposition 

N 

V>(0) = Y,Wj(0,x;a(0)) + R 
3=1 

the new perturbation ^ -J=r ^ is orthogonal to the root spaces N* of H*(a(0)) and satisfies the smallness 
condition (2.9). 

Proof. We need to solve the equation (with the solution being <r(0)) 

N 

(11.9) (v(0) -^2Wj(0, -;<7(0)), Jd Or W e (0, -,<r(0))} = for all 1 < r < 2n + 2, 

i=i 

where J is the matrix from Proposition 4.1 and Wj are the complexified versions of ip,Wj, respec- 
tively. One solves (11.9) by means of the implicit function theorem. Indeed, the derivative of the 
left-hand side of (11.9) is given by 

N N 
3=1 ' 3=1 

The second term is 0(e) where e controls the size of the initial perturbation Rq in L 1 , say. On the 
other hand, the first term is separated from zero by virtue of the convexity condition, and either the 
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separation condition (2.17) or the assumption of large relative velocities of the solitons (2.18). Indeed, 

as in Proposition 1.3.1 from [BP1] one sees that 

(11.10) 

det{(d ak W e (0,-;a(0)),Jd ar W e (0,-;a(0)))} > \\<f>(-; a,(0))||f {d a U(-; a,(0))|||) 2 > 

uniformly in the small neighborhood of a^(0) that we are allowing, see the convexity condition (2.16). 
The remaining entries of the derivative matrix, which involve inner products with Wj,Wk for j ^ k, 
are small because of either the (physical) separation condition (2.17) or the velocity condition (2.18). 
The latter ensures that we are taking scalar products of quantities that are almost orthogonal by virtue 
of the large distances of their Fourier transforms. Hence the determinant of the derivative is essentially 
bounded below by the product of the matrices with j = k, see (11.10). This proves that the derivatives 
are invertible, and since the original perturbation Rq is small, the image of the diffeomorphism given 
by the left-hand side of (11.9) contains zero, as claimed. □ 

In general, we define 
(11.11) 

id t ZW + Hia^it^Z^ = t^+^Wia^it)) + G(Z ( - n ~ 1 \a ( - n \t)) 

Z<- n \0,x) = Z (x) 

(11.12) 

s^ 1 )^^)),^^;^^ 

Here 

G{Z^ n ~ 1 \a^ n \t)) = 0(w^w^) + 0(w^w^)Z^ n ~ 1 ^ + 0(|w/( n )| p-2 |.Z'( n_1 )| 2 ) + 0{\Z^ n ~ v> \ p ) 



and 0™(i, x; a^ n \t)) once again is defined via the equation 



(»»)( 



■*j y, >~ yjj i " y \ //sj v> ' vjj ^ V! vjj < " u j 

and has the form 

Observe that by the same argument as in the case of Z^\ the perturbation Z^ is orthogonal to the 
functions ^(t,x;a^ n '(t)). We shall assume that solutions Z^ n \a^ of (11.11), (11.12) have already 
been constructed and we now proceed to estimate them. We will return to the issue of constructing 
those solutions at the end of this section. We shall now assume that aW and Z (n ~ 1 ^ satisfy (11.3) and 
(11.4) and prove the estimates (11.1) and (11.2) for er( ra+1 ) and Z^ n \ First we estimate <5-( n+1 ) in terms 
of Z^ n \ Observe that <j( n+1 ) verifies the system of ODE's described in Lemma 8.2 with the function 
G(Z( n )(t),ff( n )(t) and 0™(t,x;a( n )(t)). Therefore, 

(ii.i3) |^ n+1 )(t)| < (i + ty n {j5 2 + c\\z^\\ 2 Xs + c\\z^\\ p Xs ) 
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We now consider the Z^> equation. First, by construction Z^> is orthogonal to ^(t,x;a^(t)). 
Observe also that Z^ satisfies the equation 

idtZ^ + H{t,a (n)0 °)Z ( ^ = ^H(t,a^°°) - H(a (n \t))^Z (n ^ + t^ + ^W(a (n) {t)) + G(Z (n - 1 \a (n) (t)) 
Corollary 7.4 implies that 

H(t, a^°°) - H(a^ n \t)) = V(t, x), 

where a smooth localized potential V has the property that sup <| T |< s ||5 7 F(t, Ollz^nL 00 ^5 5 2 (l+t)~ n+2 . 
Moreover, the calculation leading to (9.6) of Lemma 9.1 shows that 

IIA^,-^, *)^ >jDj (t)^(t, -; «rW (t)) - ^*(-; o-^ 00 )!!^^ < ^(1 + t)"»+ 2 

since is an admissible path satisfying the bootstrap assumptions. Therefore, by the results of 
Corollary 5.4, taking into account smallness of the initial data Z^ n \0) = Zq, 

s 

\\Z {n) \U < HV^oIIlW + \\G{Z^\a^)\\y s + \\t^W(a^ (t))\\y s 

k=0 

(11.14) <6 2 + \\Z^f Xs + \\Z^f Xs + \\Z^\\ Xs + 

(11.15) <5 2 + \\Z^\\ 2 Xs + \\Z^\\ p Xs , 

where the inequality leading to (11.14) follows from the estimates on the nonlinear term G{Z <yU ~ 1 ^ (t),a^ (t)) 
and t^ l+1 ^W{a (n \t)) obtained in Lemma 9.3. The bound on t^ n+l ^W {a (n \t)) also uses the inequality 
(11.13). To pass to (11.15) we used the assumption that Z^ 1 ^ satisfies (11.2). In the same way one 
obtains a local in time version of equation (11.15). 

(11-16) \\Z^\\ MT) <S 2 + \\Z^\\% s(T) + \\Z^\f Xs(T) , 

which by continuity in T implies the desired estimate ||Z^ n ^||^ s < 5 2 . 

Thus, the sequence Z^ is uniformly bounded and small in the space X s while (1 + t) n d^ is 
uniformly small point wise in time. Therefore, we can choose a convergent subsequence of the paths 
a( k \t) — ► a(t) and a weekly convergent in H s (R n ) subsequence Z^ — > Z. We multiply the equation 
(11.11) by a smooth compactly supported function ((x), integrate over the entire space and pass to 
the limit using that on any compact set Z^ — > Z strongly in H s for any s' < s. In particular, since 
s > t-, Z( n ) — ► Z point wise. It will follow that Z is a solution of the equation 

(11.17) id t Z + H(a(t))Z = tW(a(t)) + G(Z(t),a(t)), 

Z(0,x)=Z (x) 

We also pass to the limit in the equation (11.12) to obtain 

(11.18) (tW(a(t)),^(t, ■; a(t))) = (G(Z(t),a(t)),^(t, •; a(t))) + (p?(t, •; a(t)), Z 
Comparing equations (11.17) and (11.18) we conclude that 

(Z(t),^(t,-;a(t)))=0 
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for all j, m. Therefore, the function %p = R + w\ + w 2 solves the original NLS and by uniqueness, say 
in L 2 , i/j is our original solution. 

To show existence of the solution Z^ n \ (j( n+1 ) of the linear system (11.11), (11.12) we first construct 
the solution on a small time interval. We note that the "system" (11.12) for <5-( n+1 ) can be resolved 

r (»)/ 

the system (11.11), (11.12) by the following caricature: 



algebraically due to the spatial separation of the paths ™ (t). Therefore, for simplicity we can replace 



idtz + = ^(*> x ) z + u{t)a{t, x ) + gii-, x )i 

Lu(t) = (z,b(t,-))+f(t) 

Here V,a,b,f are sufficiently smooth given functions and git,-) € H s uniformly in t. We eliminate 
u(t) and infer that 

id t z + - Az = V(t, x)z + {z, c(t, -))a(t, x) + F(t, x) 

with some new smooth functions c,a and an H s function F(t, •). Using the standard energy estimates 
we obtain that 

\\z(t)\\ H » < \\zo \\h° + C x t sup ||^(r)||ij- s +C 2 t sup ||F(r, -)||h*, 

T<t T<t 

where the constants C\ and C2 depend on V, a, b, f. Therefore, we can establish the existence of the 
solution on the time interval of size \C^ 1 by means of the standard contraction argument. Then we 
can repeat this argument indefinitely thus constructing a global classical solution. 

12 The linearized problem 

12.1 Estimates for matrix charge transfer models 

In this section we recall some of the estimates from Sections 7 and 8 from our companion paper [RSS]. 
First, consider the case of a system with a single matrix potential: 

*0 + (V-;%)(S)- 

with U, W real- valued and H = \ A - fi, fi > 0. We say 
iff the conditions of the following Definition 12.1 hold. 

Definition 12.1. Let A be as above with U,W real-valued and exponentially decaying. The operator 
A on Bom(A) = H 2 (R n ) x H 2 (R n ) C Ti := L 2 (R 3 ) x L 2 (R 3 ) is admissible provided 

• spec(^4) C K and spec(^4) n (— fx, fi) = {uJe I < £ < M}, for some M < 00 where ujq = and all 
ujj are distinct eigenvalues. There are no eigenvalues in spec ess (^4) = (—00, — fj] U \p, 00). 

• For 1 < t < M, L t := ker(A - uj e ) 2 = ker(A - u t ), and ker(A) C ker(A 2 ) = ker(A 3 ) =: L . 
Moreover, these spaces are finite dimensional. 
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• The ranges Ran(A - u e ) for 1 < i < M and Ran(A 2 ) are closed. 

• The spaces are spanned by exponentially decreasing functions in Ti (say with bound e~ £ °\ x \ ). 

• The points ±/x are not resonances of A. 

• All these assumptions hold as well for the adjoint A* . We denote the corresponding (generalized) 
eigenspaces by L* f . 

We will discuss these conditions in detail in the following Subsection 12.2. It is possible to establish 
some of these properties by means of "abstract" methods (for example, the exponential decay of 
elements of generalized eigenspaces via a variant of Agmon's argument, or the closedness of Ran(A— ug) 
from Fredholm's theory), whereas others can be reduced to statements concerning certain semi- linear 
elliptic operators L + ,L_, see (12.31) (for example, that the spectrum is real or that only can have 
a generalized eigenspace). In a later section we will prove for a particular model that L + , L_ have the 
required properties. One condition that we will not deal with in this paper is the absence of embedded 
eigenvalues in the essential spectrum. This property will remain an assumption. 

It is shown in [RSS], Lemma 7.2 that under these conditions there is a direct sum decomposition 

M M 



;i2.2) w = E L *+(E L J)" 



j=o 3=0 

^M 



and we denote by P s the induced projection onto ^X^j=o ^j) ■ ^ n g enerai > Ps is non-orthogonal. The 
letter "s" here stands for "scattering" (subspace). It is known that Ran(P s ) plays the role of the 
scattering states for the evolution e ltA . Indeed, the main result from Section 7 in [RSS] is that if A is 
admissible and the linear stability property 

(12.3) sup||e iM P s || 2 ^2 < oo 

t 

holds, then one has the dispersive bound 

(12.4) ¥ tA PsM\L-> +L °° < \A~Hm\l^ 



(if in addition ||V"||i < oo, then the L? norm can be removed on the left-hand side). Next, we recall 
the notion of matrix charge transfer models from Section 8 in [RSS]. 

Definition 12.2. By a matrix charge transfer model we mean a system 
(12.5) idS+(^ _? A V+X>(-"^=0 

ip\t=o = ipo, 

where Vj are distinct vectors in M 3 , and Vj are matrix potentials of the form 



Vj{hX) ~\e- ie ^Wj{x) -Uj(x) J 
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where 9j(t, x) = (\vj | 2 + a 2 )t + 2x ■ Vj + jj, aj, jj G M ; aj 7^ 0. Furthermore, we require that each 

±A - \a 2 + JJ _ w 



be admissible in the sense of Definition 12.1 and that it satisfy the linear stability condition (12.3). 

It is clear that the Hamiltonian in (2.23) is of this form. As in Lemma 3.1 above one now verifies the 
following. The Galilean transforms Q$ := Gvo are defined as in (3.11), i.e., 



fi\ = ( Qv,o(t)h \ 
h) 



where Q^ (t) = e~ i 2 ' e~ ix ' v e ito > . 
Lemma 12.3. Let a G M. and set 

A -\ W -iA + ±a 2 -U 

with real-valued U, W . Moreover, let v G M 3 , 9(t, x) = (\v | 2 + a 2 )t + 2x ■ v + 7, 7 G M., and define 

H(t\ =( l A + U (-~ *>t) -e id ^-^W(- - vt) 
[ ) : ~ \ e-W'-^Wi- - vt) -±A - U(- - vt) 

Let S{t), 5(0) = Id, denote the propagator of the system 

id t S{t) + H(t)S{t) = 0. 

Finally, let 

( e -i«(t)/2 \ 

(12.6) M(t) = M aa {t) = [ Q e 4 )/2 J 
where u>(t) = a 2 t + 7. TTien 

(12.7) S(t) = g i r,o(i)" 1 ^W~ 1 e' M A^(0)^,o(0). 
Proof. One has 

(12.8) idtMWGifWSit) = ( ^ _ 1 ^ ) M(t)g 9 (t)S(t) + MWWijWSit) - M{t)Qz{t)H{t)S{t). 
Let p(t, x) = i|u | 2 + 2x • v. One now checks the following properties by differentiation: 

M(t)iG,(t) = _^W\ 2 + v-P _^ |2 ° + ^ ^M{t)G v {t) 

M(t)9,(t)H(t) = ( e J { tt") w li A - (7 J 

(12.9) ..p^.^AlMftM. 
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The right-hand side of (12.9) arises as follows. First, the Galilei transform introduces a factor of e 



which needs to be commuted with ^A. Since 



= \\v\ 2 e- ix -*f - iv-v(fe- ix - ff ) + \e~ ix -*M 

= (\\v\ 2 + v-pj(fe-^) +\e~™- *Af, 

one obtains the final term on the right-hand side of (12.9). It remains to check the terms involving 
the potentials (for simplicity 9(- — tv) = 9(t, ■ — vt)): 



M{t)g,{t) 



U(- - Vt) - e iO(t,-vt) W ^ _ ^ \ /f x 

e -iO{t,-vt) W ^ _ ^ _ ^ J ^ 



e 



iw(t)/2 o \ f^(t)U{- - vt)h - ^{t)e^-^W{- - vt)f 2 \ 



e^ W/2 J \Uty e{ -^W{- - vt)h ~ Qv(t)U(- ~ vt)h) 



We i(vH+2x-v) e i(u,-e) te (t)(e-^(*)/ 2 /!) - U 9v(t)e^( t )/ 2 f 2 ) 



-UJ- 



p)W -U J V e*"<*>/ 2 J \ Q ,(t)J 



as claimed. In view of our definitions, 9 — p—uj = 0. Since uj = a 2 , the lemma follows by inserting (12.9) 
into (12.8). □ 

In order to prove our main dispersive estimates for such matrix charge transfer problems we need 
to formulate a condition which ensures that the initial condition belongs to the stable subspace. To do 
so, let P s (Hj) and P b (Hj) be the projectors induced by the decomposition (12.2) for the operator Hj. 
Abusing terminology somewhat, we refer to Ran(Pb(Hj)) as the bound states of Hj. 

Definition 12.4. Let U(t)tpo = ip(t, ■) be the solution of (12.5). We say that ipo is a scattering state 
relative to Hj if 

\\P b (Hj,t)U(t)ifo\\ L 2 ^0 as t^+oo. 

Here 

(12.10) Pb(Hj,t) := g ilj (t)~ 1 M :i (t)- 1 P b (H J )M j m ; ; j (t) 

with Mj(t) = M ajnj (t) as in (12.6). 

The formula (12.10) is of course motivated by (12.7). Clearly, P b (Hj,t) is the projection onto 
the bound states of Hj that have been translated to the position of the matrix potential Vj(- — tiJj). 
Equivalently, one can think of it as translating the solution of (12.5) from that position to the origin, 
projecting onto the bound states of Hj, and then translating back. 

We now formulate our decay estimate for matrix charge transfer models, see Theorem 8.6 in [RSS]. 
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Theorem 12.5. Consider the matrix charge transfer model as in Definition 12.2. Let U(t) denote the 
propagator of the equation (12.5). Then for any initial data ipo G L 1 n L 2 , which is a scattering state 
relative to each Hj in the sense of Definition 12.4, one has the decay estimates 

(12.11) \\U{t)^\\ L oo <{t)-HMWnL*- 

For technical reasons, we need the estimate (12.11) for perturbed matrix charge transfer equations, 
as described in the following corollary. This is discussed in Remark 8.6 in [RSS]. 

Corollary 12.6. Let ip be a solution of the equation 

(12.12) _° lA ^+J2v j (--v j t)i;+V (t,x)^=0 

V 2 / j=1 

ip\t=o = ipo, 

where everything is the same as in Definition 12.2 up to the perturbation Vo(t, x) which satisfies 

sup ||Vb(t, -)\Wr\L^ < £■ 
t 

Let U(t) denote the propagator of the equation (12.12). Then for any initial data ipo € L 1 n L? , which 
is a scattering state relative to each Hj in the sense of Definition 12.4 (with U(t) replaced by U(t)), 
one has the decay estimates 

(12-13) \\U(Wo\\ L ~<(t)-l\\j \\ LlnL2 
provided e is sufficiently small. 

The corresponding inhomogeneous bound is stated in Section 12. 



12.2 The spectral properties I: general arguments 

In order for the linear estimates to apply, we need to impose the conditions in Definition 12.1 as well 
as the linear stability condition (12.3) on the operators from (2.24). The admissibility conditions of 
Definition 12.1 were motivated to a large extent by Buslaev and Perelman [BP1], who built on earlier 
work of Weinstein [Wei]. We now analyze these conditions in detail. As before, 

< 12 - 14 ) A --={ W -H-U )= B + V 

where U, W are real- valued, H = \ A — fi with fi > 0, and V is the matrix potential consisting of U, W. 
In this subsection we deal with those properties that can be dealt with by means of general arguments, 
that make no use of any special structure of the operator. 

Lemma 12.7. Let the matrix potential V be bounded and go to zero at infinity. Then (A — z) 1 
is a meromorphic function in Q := C\ (— oo, — fi] U \p,oo). The poles are eigenvalues of A of finite 
multiplicity and Ran(yl — z) is closed for all z G Q,. Finally, the complement of Q agrees with the 
essential spectrum of A, i.e., spec ess (yl) = (— oo, — ji\ U [/x, oo). 
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Proof. Suppose that z£(l. Then B — z is invertible, and A — z = (l + V(B — z) ^ (B — z). Since 

V(B — z)~ x is analytic and compact in that region of z's, the analytic Fredholm theorem implies that 
1 + V(B — z)' 1 is invertible for all but a discrete set of z's in Q. Furthermore, the poles are precisely 
eigenvalues of A of finite multiplicity. It is also a general property that the ranges Ran(l + V (B — z)~ v ) 
are closed. Indeed, if K is any compact operator on a Banach space, then it is well-known and also 
easy to see that Ran(J — K) is closed. Since B — z has a bounded inverse for all zeSl, this implies that 



Ran(^4 — z) is closed, as claimed. Conjugating by the matrix P 
(12.15) 



A ■= P- X AP = i 



H + Vi 
-H-V 2 



B + V, B = i 





i 


(! 


—i 





H 


-H 






leads to the Hamiltonians 



V 



Vi 
-V 2 



where V\ = U + W and Vi = U — W. The system (12.15) corresponds to writing a vector in terms 
of real and imaginary parts, whereas (12.28) corresponds to working with the solution itself and its 

i 



conjugate. By means of the matrix J 



B 



H 
H 



-i 
J, A 



one can also write 



H + Vi 
H + V 2 



J. 



Since B* 



B it follows that spec(i?) C I 

(H 2 



(B-z)- 1 = 



(12.16) 
(12.17) 



(B + z) 
{H 2 - 



One checks that for Rz ^ 

z 2 )- 1 







(H 2 - z 2 



-l 



-z 2 )- 1 
(H 2 -z 2 )- 1 



{A-z)- 1 = (B-z)- L -(B-zy'W! 



-i 



(B + z) 
1 + W 2 J(B - z)~ l Wi 



-i 



W 2 J{B - z) 



-i 



where W\ and W 2 are the following matrix potentials that go to zero at infinity: 







\V 2 



w 2 



|Vi|2sign(Vi; 






|V r 2 |2sign(V2) 



The inverse of the operator in brackets exists if z = it with t large, for example. Moreover, by the 
assumed decay of the potential the entire operator that is being subtracted from the right-hand side 
is compact in that case. One is therefore in a position to apply Weyl's criterion, see Theorem XIII. 14 
in [RS4], whence 



(12.18) spec ess (A) = spec css (I) = (-oo, -//] U [//, oo). 

The identity (12.17) goes back to Grillakis [Gr]. 



□ 



Next, we need to locate possible eigenvalues of A or equivalently, A. This will not be done on the 
same general level, but require analysis of L + ,L_ from (12.31). But we first discuss another general 
property of the matrix operator A. 
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Lemma 12.8. Let A be as in (12.14) with U,W continuous and W exponentially decaying, whereas U 
is only required to tend to zero. If f G ker(^4 — E) k for some —p < E < p and some positive integer k, 
then f decays exponentially. 

Proof. We want to emphasize that the following result is "abstract" and does not rely on any special 
structure of the matrix potential or on any properties of L + or L_. We will use a variant of Agmon's 
argument [Ag]. More precisely, suppose that for some —p < E < p, there are ^1,^2 & -£f 2 (M n ) so that 

(A - p + U)ip! - Wi> 2 = Efa 
(12.19) Wtpi + (- A + p - U)tp2 = E^ 2 . 

As usual, U,W are real-valued and exponentially decaying, fi > 0. Suppose |W(x)| < e -6 l x L Then 
define the Agmon metrics 



p%{x) = inf L± ( 7 ) 

7:0— s 



(12.20) L± g ( 7 ) = f 1 mm(y/( f i±E-U( 1 (t))) + ,b/2) 

JO 



dt 



where j(t) is a C 1 -curve with t G [0, 1], and the infimum is to be taken over such curves that connect 
0, x. These functions satisfy 



(12.21) |Vp|(x)| < y/(ii±E-U{x)) + . 

Moreover, one has p%{x) < b\x\/2 by construction. Now fix some small e > and set w ± (x) := 
e 2(l-e)p B (x)_ Q ur g oa j j g to g^Q-^r that 

(12.22) / L+^lM^l 2 +u-(x)\i/> 2 (x 



dx < 00. 

Not only does this exponential decay in the mean suffice for our applications (cf. Section 7 in [RSS]), 
but it can also be improved to point wise decay using regularity estimates for ipi,ip2. We do not 
elaborate on this, see for example [Ag] and Hislop, Sigal [HiSig]. 
Fix R arbitrary and large. For technical reasons, we set 



p% R (x) := min (2(1 - e)p±(x), i?) , w±(x) := 



Notice that (12.21) remains valid in this case, and also that p%{x) < min(6|x|/2, R). Furthermore, by 
choice of E there is a smooth functions cj) that is equal to one for large x so that 

supp(^) C {p + E - U > 0} n {p - E - U > 0}. 

It will therefore suffice to prove the following modified form of (12.22): 

(12.23) sup / ru+(x)|-(/;i(x)| 2 +uj R {x)\ip 2 {x)\ 2 \4> 2 {x)dx < 00. 

R J L - 1 
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All constants in the following argument will be independent of R. By construction, there is 5 > such 
that 

(12.24) 5 J uj+(x)^ 1 (x)\ 2 (f) 2 (x)dx < J uj+(x)(fi + E -U(x))^ l (x)\ 2 <p 2 (x)dx 
= J u^(x)(A^i -W^ 2 )(x)Mx)^ 2 (x)dx 

(12.25) = - j V(u+(x)0 2 (x))VMx)Mx) dx - J 4(i)^ 2 (i)|V^(i)| 2 dx 

(12.26) - J u^{x)W{x)^ l {x)i) 2 {x)<j) 2 {x)dx. 

As far as the final term (12.26) is concerned, notice that sup^ \u^(x)(/) 2 (x)W(x)\ < 1 by construc- 
tion, whence | (12.26) [ < HV'ilbllV^lh- Furthermore, by (12.21) and Cauchy-Schwarz, the first integral 
in (12.25) satisfies 



(12.27) <2(l-e)(y oj+(x)(h + E - U(x))^ 2 (x)\Mx)\ 2 dx} 1 (J u+{x)ct){x) 2 |V^i(a;)| 2 dx} 1 

+2(y lj+{x)4> 2 {x)\v^{x)\ 2 d X y (y c^iv^i 2 i^i 2 d X y 

Since the first integral in (12.27) is the same as that in (12.24), inserting (12.27) into (12.25) yields 
after some simple manipulations 

e J Lo+(x)(fi + E -U(x))^ 1 (x)\ 2 <p 2 (x)dx < e' 1 J uj+(x)\V<f)(x)\ 2 |V>i(x)| 2 dx 

— j u!^(x)(j)(x) 2 W(x)ip 2 (x) ipi(x) dx. 



Since Vcfi has compact support, and by our previous considerations involving uj^W, the entire right- 
hand side is bounded independently of R, and thus also (12.24). A symmetric argument applies to 
the integral with ip 2 , and (12.23), (12.22) hold. This method also shows that functions belonging to 
generalized eigenspaces decay exponentially. Indeed, suppose (A — E)g = and (A — E)f = g. Then 

(A - ft + U)h ~ Wf 2 = Ef 1+9l 
Wh + (- A + n - U)f 2 = Ef 2 +g 2 

with g±,g 2 exponentially decaying. Decreasing the value of b in (12.20) if necessary allows one to use 
the same argument as before to prove (12.22) for /. By induction, one then deals with all values of k 
as in the statement of the lemma. □ 

12.3 The spectral properties II: reduction to L + ,L 

We now need to specialize A from (12.14) to the form (2.24), i.e., 

n o - A = (\ A " £ + /3(0 2 ) + W 2 )0 2 /W)<A 2 \ 

1 ] \ -/?'(<w -*a + £ - am - /n w ) • 
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As shown in Section 2, these are the stationary Hamiltonians derived from the linearization of NLS, 
see (2.24) and Lemma 12.3. Let a > and be a nonzero solution of 

(12.29) l A 0-^ + /3(0 2 )0 = O, 

i.e., 4> = (/)(■; a) for all a G (ao — cq, «o + cq). Moreover, <p is smooth in both variables and 

(12.30) ||0a0||tfi(R») + Wdl^WuH^) < oo. 

Finally, we require that 4> is sufficiently rapidly decaying. A particular case would be a which is 
positive and radially symmetric. Such a solution is known to exist and to be unique if (3{u) = \u\ a 
provided < a < an d is referred to as the "ground state". It decays exponentially. However, in 
order to keep this section as general as possible, we do not require 4> to be the ground state. Let 

(12.31) L := ~A + ^ - /3(<A 2 ), L + := ~A + ^ - /3(0 2 ) - 2(3' ' 
with domains Dom(L + ) = Dom(L_) = H 2 (W l ) so that 



(12.32) 




withDom(Ji) = H 2 (M. n ) x H 2 (R n ). Here A is obtained by conjugating A with the matrix P, see (12.15). 
For simplicity, however, we no longer distinguish between A and A, i.e., we set A = A. The spectrum 

2 

of L± on [fi, oo) is purely absolutely continuous, and below [i = ^- > there are at most a finite 
number of eigenvalues of finite multiplicity (by Birman-Schwinger, the assumed decay of <f> as well as 
0(0) = 0). Clearly, 

(12.33) L_0 = O, L + (dj<f>) = 0, 1 < j <n, L+(d a <l>) = -a<f>, 

where the final property is formal. We now collect some crucial properties discovered by M. Weinstein. 



Definition 12.9. 



Needed properties of the scalar elliptic operators L + and L ; 



Let(fr(-;a) be as above, in particular assume that (12.30) holds. The kernels have the following explicit 
form: 

ker(L_) = span{</>} and ker(L + ) = span{<9j0 | 1 < j < n}. 

The operator L + has a single negative eigenvalue E\ with a unique ground state ip > 0, whereas L_ is 
nonnegative, with as an isolated eigenvalue. Furthermore, the convexity condition (d a (p(-; a), a)) > 
holds, see (1.17). 



In the following section, these conditions will be verified for a particular choice of nonlinearity. 
Note that according to this definition, (j) is the ground state of the linear operator L_ and as such is 
positive. In case of 4>(-;a) being the ground state of the nonlinear problem (12.29), these properties 
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have been shown to hold by Weinstein [Wei] and [We2] in case of power nonlinearities, i.e., (3(u) = \u\ a , 

°<^< uk- 

The purpose of this subsection is to reduce some of the admissibility conditions from Definition 12.1 
to the properties of L+, L_ from Definition 12.9. Recall from the previous subsection that several other 
properties hold in greater generality. We now collect those properties that follow from Definition 12.9 
into a single proposition. 

Proposition 12.10. Impose the spectral assumption on L + and L_ from Definition 12.9. Then 

• spec(A) C R, the only eigenvalue that admits a generalized eigenspace is 0, and R&n(A 2 ) is 
closed. 

• the linear stability condition holds. 

• equality holds in the relation concerning the root spaces (12.39) and (12.40). In particular, one 
has ker(^4 2 ) = ker(^4 3 ) and ker((^4*) 2 ) = ker((A*) 3 ) and the rootspace J\T(A*) has dimension 



The proof of this proposition is split into several lemmas below. 

Lemma 12.11. Impose the spectral assumption on L + and L_ from Definition 12.9. Then spec(yl) C 
1R, the only eigenvalue that admits a generalized eigenspace is 0, and Ran(^4 2 ) is closed. 

Proof. Consider 



with domain H (W 1 ) = W 4,2 (W l ). Following [BP1], we first show that any eigenvalue of T, and 
therefore also of spec(yl 2 ) is real, and then under the assumption (1.17), that it is nonnegative. Because 
of (12.18), the latter then implies that spec(A) is real, as required in Definition 12.1. Clearly, T(f> = 0. 
Let Y> ^ span{(/>}, Tip = Eip. Let ip = ipi + ccf), tpi _L <j). Then 



so that L^ipi ^ is an eigenfunction of the symmetric operator L£L + Ll (with domain H 4 (W 1 )), and 
thus E is real. Hence any eigenvalue of A can only be real or purely imaginary. Since (j) _L ker(L + ) by 
our assumption concerning L + , the function 



2n + 2. 



(12.34) 
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in view of (12.33) and (1.17). Now suppose that A 2 has a negative eigenvalue. Then by the preceding, 

i i 

so does T, and therefore also LiL+Li. More precisely, the argument from before implies that there 
is x £ ker(L_)- L , x / 0, so that 

(l!l + l!x, x ) = (l + v>,V)<o 

i 

with ip = L'ix- Let Pi denote the projection onto the orthogonal complement of ker(L_) = span(^). 
By the Rayleigh principle this implies that the self-adjoint operator P^L + P^ has a negative eigenvalue, 
say £3 < 0. Thus L + tp = E^ip + op for some tp _L (p. If c = 0, then £3 = E\ so that ip > as the 
ground state of L + . But then ((p,ip) > 0, which is impossible. So c / 0, and one therefore obtains 

(L+ - £ 3 ) _ V = -J> g(E 3 ) = 0. 

But this contradicts (12.35) by strict monotonicity of g. Thus A 2 does not have any negative eigen- 
values, which implies that A does not have imaginary eigenvalues. Hence all eigenvalues of A are real, 
as desired. 

We now turn to generalized eigenspaces. Suppose Atp = Etp + x, where E ^ 0, (A — E)x = and 
X 7^ 0. This is equivalent to saying that A has a generalized eigenspace at E. Then G Dom(^4 2 ), 
and moreover 

(A 2 - E 2 ) X = 0, (^ 2 - E 2 )^ = (A - E) X + 2L X = 2E X , 

so that A 2 would have a generalized eigenspace at E, and therefore also T. Hence, suppose Tip = Etp, 
with E + 0, ip ^ 0. If (T — £)x = with c / 0, then 

{LlL+Ll - E)Ll X i = cLlip! + 0, {LlL + Ll - E) 2 h\x\ = cLl(L + L_ - E)ip = 

1 

where ipi,Xi denote the projections of tp, x onto the orthogonal complement of <p. But Ltip\ 7^ since 

1 1. 

E / and thus E would have to be a generalized eigenvalue of LiL + Li, which is impossible. So only 

1 1 

E = can have a generalized eigenspace. Here we used the property that LiL + Li is self-adjoint on 
its domain Lf 4 (IR n ). While symmetry is obvious, self-adjointness on H A (W l ) requires a bit more care. 

Suppose {LlL+Llf,g) = (f, h) for all / G H 4 (R n ), and some fixed g,he L 2 (R n ). Taking / G ker(L_) 

1 

shows that Pi/t = h, i.e., that h € (ker(Li))- 1 . By the Predholm alternative applied to the self-adjoint 

11 1 
operator Li, one can write h = L 2 _h\ with some h\ G Dom(Li) = H 1 (W l ). Note that h\ is defined 

1 

only up to an element in ker(Li), i.e., hi + ccj> has the same property for any constant c. Thus 

(LlL + Llf, g) = (/, Llih + af>)) = {L*f, h ± + af>) 

1 

for all / G Lf 4 (R n ). Equivalently, setting /1 = Li/, one has 

{LlL + h,g) = {hM + c<j>). 
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Note that the class of /i are all functions in H 3 (M. n ) with fi _L (ft. We now want to remove the latter 
restriction, which can be achieved by a suitable choice of c. Indeed, in order to achieve 

{LlL + (fx + X(f)),g) = (h + ^M + ccj)) 

for all fi G H 3 (M. n ), fi _L (ft, A G C one chooses c such that 

i 

(L?_L + <ft,g) = ((f), hi + c(ft), 
which can be done since ((ft, (ft) > 0. Renaming h\ + ccft into h±, one thus arrives at 

(12.36) (LlL+f^g) = (hM) for all f x G H 3 (R n ). 
i 

Recall that /i = L1h\. One can now continue this procedure. Indeed, since (12.36) implies that 
h\ _L ker(L + ), one can write h\ = L + /i 2 = L + {h,2 + ]Cj=i c jV':/)> where ker(L + ) = span{V>j}j =1 and 
/i 2 € -fT 3 (M n ) (in fact, iftj = dj(ft by our assumption). As before, the constants {cj} are chosen in such 
a way that 

A; fc fc 

for all Xj. This can be done because of the invertibility of the Gram matrix of {ipj}j =l . Hence 

(Llf 2 ,g) = (f 2 ,h 2 ) for all /^ff 1 ^"). 

i ii 
Moreover, h = L'tL + h 2 with h 2 G H 3 (R n ). By the self-adjointness of L 2 _ this implies that hi = Ltg. 

1 JL 

It follows that g G F 4 (R n ) and h = L 2 _L+L 2 _g as desired. 

Finally, we show that Ran(^4 2 ) is closed. By (12.34) it suffices to show that the ranges of both 
T = L + L^ and T* = L_L + are closed with domain H A (R n ). We will first verify that these operators 
are closed on this domain. Indeed, they each can be written in the form A 2 + F\ A + AF 2 + F3. Since 
for M large 

||(A 2 + M + Fi A + AF 2 + F 3 )/|| 2 > ||(A 2 + M)/|| 2 - C\\f\\ w ^ > i||(A 2 + M)f\\ 2 > \\f\\ w ^, 

one concludes that T + M, T* + M are closed, and therefore also T,T*. Next, let P_ and P+ be the 
projections onto ker(L_) and ker(L + ), respectively. Then L_L + = L-P^L + P+ . Now 

(12.37) P^L + f = L+f - \\(ft\\- 2 (L + f,(ft)(ft = L+(/ - IMpV, L+<^) , 

where we have written = G ker(L + )^ = Ran(P^r) by virtue of the fact that 

(ft G Ran(L + ) = Ran(L + ) = kcr(L + )" L = spanlc^} -1 . 
The last equality here is Weinstein's characterization, more precisely, our assumption on L + . Define 
4>:=P^i, Qf:=f-\\(ft\\- 2 (f,L + (ft}4>, Qf:=f-\\(ft\\- 2 {f\PiL + (ft)^ 
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so that (12.37) gives 



P-L+P^ = L + QP^ = L + P^Q => L_L + = L_L+P^Q 



In particular, 



||T7|| 2 = ||L_P±L + P|/|| 2 > Cl ||P±L+i^/|| 2 
(12.38) = Cl ||L + QP|/|| 2 =ci||L + P|Q/|| 2 > c lC2 \\P^Qf\\ 2 , 

where the existence of ci,c 2 > follows from the self-adjointness of L_,L + . Hence, if T*f n = 
T*P^Qf n -> /i in L 2 , then by (12.38) and linearity P+Qf n -> £/ in L 2 . Since T* was shown to 
be closed, it follows that h = T*g, and Ran(T*) is closed. A similar argument shows that Ran(T) is 
closed. □ 



Next, we derive the linear stability assumption as well as the structure of the generalized eigenspaces 
of A and A* from our spectral assumptions on L_, L + . From the spectral assumptions in Definition 12.9 
as well as 

L_0 = 0, L^(xj4>) = —dj(j), L + (dj(p) = 0, L + {d a 4>) = —a<\> 

it follows that 

ker(A) = span{("), (ff) ■ 1 < J < »} 
ker (j f) = spa„{(*). Q :!<,<„} 

(12.40) =- UM(^)')3s P a„{(*), (^), (f), ( 8 °J :l<i <„}=:*.. 

One of our goals is to show that equality holds in the last two relations. This is the same as the 
structure statement made in Proposition 4.1, but one needs to apply the matrix P = ^ J * . ^ to 
pass between these two representations. 

Now suppose that zdtV + = 0. This can be written as dtip + JMip = where J = 

and M = ^ ^° ^ . Therefore, 

^(tf, M$) = 23?(d t ^, Mip) = -25ft( JM^, Mip) = 
by anti-selfadjointness of J. In other words, 

Q$) := (L + ^i,Vi) + (^2,V'2) 



52 



is constant in time if tf(t) = e itA ip(0) (here $ = (J)). Although the previous calculation basically 
required classical solutions, it is clear that its natural setting are ii -1 (Resolutions. In that case one 
needs to interpret the form via 

(L + Vi,Vi) = |l|V^i||l + ^H^illl - </5f(0 2 )^i,^i) 

(12.41) (L-M2) = i||V^ 2 ||2 + ^||^ 2 ||2-((/3(0 2 ) + 2/3 / (0 2 )^)^,^2). 

In what follows, we will tacitly make this interpretation whenever it is needed. The following lemmas 
are due to Weinstein [Wei]. 

Lemma 12.12. Impose the spectral assumptions on L + ,L- from Definition 12.9. Then (L + f,f) > 
for all f eH^W 1 ), f ± <j>. 

This is a special case of Lemma E.l in [Wei], and we refer the reader to that paper for the proof. 

Lemma 12.13. Impose the spectral assumptions on L + ,L- from Definition 12.9. Then there exist 
constants c = c(a,(3) > such that for all ip £ H 1 (M n ), 

1. (£-^2,^2) > c\\ip 2 \\l if i>2 -L d a 4>, ^ 2 -L djcj) 

2. (L+^i^i) > c|| Villi if ipi -L 4>, ipi -L Xj(f>. 

The constant c(q, (3) can be taken to be uniform in a in the following sense: If «o satisfies Defini- 
tion 12.9, then there exists 5 > so that 1. and 2. above hold for all \a— ao\ < 5 with c(a, [3) > ^c(ao, (3). 

Proof. Consider the minimization problems 

(12.42) ini^L-f, f) subject to constraints ||/|| 2 = 1, / JL d a (f), f J_ dj<f> 

(12.43) inf (L + g,g) subject to constraints \\g\\2 = 1, g -L <j>, g -L Xj<fi. 

As usual, one would like to establish the existence of minimizers by means of passing to weak limits in 
minimizing sequences. While such sequences are bounded in H 1 (M n ), this is not enough to guarantee 
strong convergence in L 2 (R n ) because some (or all) of the L 2 -mass might escape to infinity. Using the 
fact that the quadratic forms in question are perturbations of 5IIV/H 2 . + ^rll/Hl by a potential that 
decays at infinity, one can easily exclude that all the L 2 -mass escapes to infinity. One then proceeds 
to show that the remaining piece of the limit, normalized to have L 2 -norm one, is a minimizer. This, 
however, is a simple consequence of the non-negativity of L_ and L + , the latter under the constraint 
/ _L (f), see Lemma 12.12 above. This argument is presented in all details in [Wei], page 478 for the case 
of power nonlinearities. But the same argument also applies to the general nonlinearities considered 
here, and we do not write it out. 

Assume therefore that /o is a minimizer of (12.42) with H/0II2 = 1, /o -L d a (f), /o _L dj4>. Then 

n 

(12.44) L_/ = A /o + c cU + £ Cjdrf 
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for some Lagrange multipliers \$,co, . . . ,c n . Clearly, Ao agrees with the minimum sought, and therefore 
it suffices to show that Ao > 0. If Ao = 0, then taking the scalar product of (12.44) with (j> implies 
that Co = (using that (d a cj),(j)) > 0). Taking scalar products with Xk4> shows that also Ck = for 
1 < k < n. Thus L_/o = 0, which would imply that /o = j(p f° r some 7 / 0. However, this is 
impossible because of fo _L d a (p. 

Proceeding in the same manner for L + , one arrives at the Euler-Lagrange equation 

n 

L+go = Xogo + c ^ + ^ CjXjcf). 

As before, Ao is the minimum on the left-hand side of (12.43) and thus Ao > by Lemma 12.12. If 
Ao = 0, then taking scalar products with dk4> leads to c& = for all 1 < k < n. Hence L + go = co<t> 
which implies that 

n 

go = d a 4> + y~] b e de<p. 

a f— ' 

i=\ 

Taking scalar products of this line with (j) and Xj(fr shows that Co = and be = for all 1 < I < n, 
respectively. But then go = which is impossible. 

Since the constants c(a, [3) > were obtained by contradiction, one has no control on their depen- 
dence on a. However, let \a — ao\ < 5 be as in Definition 12.9. Suppose ||/||2 = 1 satisfies / ± cj)(-, a), 
f J_ Xj(f)(-,a). Then there is 

h £ sp&n^(f)(-,a), (f)(-,ao), Xj(/)(-,a), Xj<j)(-,ao) : 1 < j < n| 

so that f + h A. (/)(■, oq) and / + h ± Xj(f>(-, ao). Moreover, since ||9 Q 0||^i( K n) + H^^H^i^n) < 00 
one can take ||/i||_f/i(Rn) as small as desired provided 5 is chosen small enough. One can therefore use 
inequality 2. from this lemma at ao for / + h to obtain a similar bound for / at a. □ 

The following corollary proves the crucial linear stability assumption contingent upon the spectral 
assumptions on L + , L_ from above (and thus, in particular, contingent upon the convexity condition). 
Strictly speaking, the following corollary gives a stronger statement than (12.3), since the range of P s 
is potentially smaller than needed for the stability to hold. 

Corollary 12.14. Impose the spectral assumptions on L + ,L- from Definition 12.9. Then there exist 
constants C = C(a, (3) < 00 so that for all 1^0 £ H 1 (W l ) 

(12.45) \\e ltA ipo\\m(R") < C||^o||ffi(R») provided tp G . 

Here M.* is the A* -invariant subspace from (12.40) . Moreover, the same bound holds for H s {W l )-norms 
for any real s with s-dependent constants (and thus in particular for L 2 (M. n )). Analogous statements 
hold for e ttA * . Finally, the constants C(a,/3) can be taken to be uniform in a in the following sense: 
If ao satisfies Definition 12.9, then there exists 5 > so that (12.45) holds for all \a — ao\ < 5 with 
C(a,P) < 2C(a ,(3). 
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Proof. Let ijj(t) = e ltA tpQ. By Lemma 12.13 one has 

Q(jo) = Q(m)>c\\m\\ 2 L2{Rn) 

provided that Vo £ -M*- Since clearly Q{i[>o) < CHV'oll^i^n) one concludes that (12.45) holds with 
L 2 (R n ) on the left-hand side. In order to pass to H 1 (R n ) write 

(L.f,f) = (l-e)(L.f,f) + £ -\\Vf\\l + ef\\f\\ 2 2 -e[ ^ 2 (x))\f{x)\ 2 dx 
(12-46) > |||V/||1 + c(l - e )||/||| + e{a 2 /2 - WU)\\f\\l 

where the constant c in (12.46) is the one from Lemma 12.13. Taking e small enough, one sees that 
the third term can be absorbed into the second. Thus the entire right-hand side of (12.46) admits 
the lower bound fll/H^-i^n)- The same argument applies to L + , and (12.45) follows. The uniformity 
statement concerning the constants C(a,(3) is an immediate consequence of the analogous statement 
in Lemma 12.13. To obtain (12.45) for all H s spaces note first that 

(12.47) C7 1 ||^||^ (R „) < \\(A + iM)^\\ 2 <C e \m H ,e {Rn) 

for all integers I and sufficiently large M = M{£). Indeed, to check the lower bound for I = 1 one can 

-l 



use (A + iM)- L = (B + iM)- 1 1 + V{B + iM) 



The inverse of the operator in brackets exists 



provided M is large and it is a bounded operator on L 2 {W l ). Taking powers of this relation allows 
one to deal with all £ > 1 (in our case V is C°° which is needed here). Since A4* is ^-invariant and 
therefore M.* is A-invariant, inserting (12.47) into (12.45) allows one to pass to all odd integers s. 
The case of general s then follows by interpolation. Finally, since all arguments in this section apply 
equally well to A* as A, the corollary follows. □ 

This corollary has an important implication concerning the structure of the root spaces as required 
in Proposition 4.1. 

Corollary 12.15. Impose the spectral assumptions on L + ,L^ from Definition 12.9. Then equality 
holds in the relation concerning the root spaces (12.39) and (12.40). In particular, one has ker(A 2 ) = 
ker(A 3 ) and ker((A*) 2 ) = ker((^4*) 3 ). 

Proof. Suppose dim(J\f (A)) > 2n + 2. Then there exists ipo G AT(A) such that ipo G M.*. This is 
because a system of 2n + 2 equations in 2n + 3 variables always has a nonzero solution. Since d a <fi JL 4> 
and dj4> JL Xj<fi, one checks that ^ ker(A). Therefore, Vo £ ker(A fc ) \ ker(yl fc_1 ) for some k > 2. 
Expanding e %tA into a series implies that He^^olb > ct fc_1 for some constant c > 0, which contradicts 
Corollary 12.14. Therefore, dim(A^(^4)) < 2n + 2. Since moreover <p > and (d a 4>, (f>) > imply that 
the 2n + 2 vectors on the right-hand side of (12.39) are linearly independent, equality must hold as 
claimed. Analogously for (12.40). □ 
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13 Uniqueness of the negative eigenvalue of L + 



In this section we show that the assumption that the linear operator L + has a single simple negative 
eigenvalue of Definition 12.9 is automatically satisfied for a large class of nonlinearities. 
We consider the ground state of the problem 

(13.1) _1 A 0-/?(H 2 )0 + ^ = O 

constructed from the constrained minimization problem for the following functional: 

(13.2) J[u] = { [ |Vu| 2 : W[u] = [ G{u) = l}, 

where G(s) = Jq (/3(s 2 )s — ^s). If w is a minimum it solves the equation 



(13.3) -^Aw - X(p{w 2 )w - y w = 

where the Lagrange multiplier A is determined from the condition that W[w] = 1. We can then find a 
ground state via rescaling 

(13.4) <p{x) = w{\^x) 

for an appropriate choice of A. In this section we prove the following: 

Theorem 13.1. Let <p be the ground state of (13.1) constructed in the variational problem (13.2) with 
the nonlinearity (3 verifying the condition that \(3{s 2 ) < C|s| p_1 with p G (1, 1 + 7^5) an d that there 
exists sq such that C(sq) > 0. Moreover, assume that 



P{s 2 )s 2 > 2 f /3{r 2 )TdT 
Jo 



for all s > 0. Then the operator L + = — — [3{<j) 2 ) — P'(4> 2 )<p 2 + ^ has a unique simple negative 
eigenvalue. 

Proof. Let us consider a smooth family of functions w z satisfying the constraint W [w z ] = 1 for all z and 
passing through a minimizer w, i.e., wq = w. Observe that since w is a minimizer of the constrained 
variational problem we have that 

d d 2 



Denote 



Then 



W = —W z \ z=0 , W = —tW z \ z=0 

az az z 

d f 
= -JK]| 2=0 = 2 / VwVw 

dZ Jjgn 
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and 



d 2 f f 

(13.5) < — J[w z ]\ z=0 = 2 / |Vu>| 2 + VwViu = -2 / Aww + Awiu 

Moreover, since the family w z verifies the constraint W[u;J = 1 we have 



Computing the derivatives we obtain 

„2 



(13.6) / {I3{w 2 )w - — w)w = 0, 

/• 2 2 

(13.7) / ((3(w 2 ) + 2(3'(w 2 )w 2 - — )\w\ 2 + {f3{w 2 )w - — w)w = 
it™ 2 2 

Now observe that the Euler-Lagrange equation (13.3) implies that 

— / Aww = 2A / {[3{w 2 )w ^rw)w 

Jw 1 Jw 1 2 

On the other hand from (13.7) we have that 

r 2 r 2 

/ {f3{w 2 )w-^-w)w = - (l3(w 2 ) + 2(3'{w 2 )w 2 - ^-)\w\ 2 
Jm™ 2 j]Rn 2 

Substituting this into (13.5) we obtain that 

r 2 

(13.8) -2/ (Aw + 2\(I3{w 2 ) + 20(w 2 )w 2 -— )w)w>0 

Recall the definition of the operator L + associated with a ground state 4>: 



1 (T 2 

L + = --A-/^ 2 )-2/n<W + T 

We know from (13.4) that </>(x) = w(A~2x). Thus rescaling (13.8) leads to the inequality 
(13.9) A? +1 (L+u>(A-i) , w(\-^)) > 

Furthermore, (13.6) implies that the function u>(A~2.) can be chosen arbitrarily from the subspace 
orthogonal to the function 

* = (/3(0 2 )0-^) 

Let n denote the orthogonal projection on the above subspace. Then (13.9) implies 4 that the operator 
nL + n > 0. Now let fi and u be correspondingly a negative eigenvalue of L + and the associated 
eigenfunction. Decompose 5 

u = v + * 



*It is easy to see that the condition that ensures that A > is that 

I3{s 2 )s 2 > 2 f !3{r 2 )Tdr 
Jo 

3 One can always achieve this decomposition by normalizing u 
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where Uv = v. Then 

L + u = L+v + L+V = n{v + 

Projecting we obtain that 

(ILL+II - n)v = -ILL+tf 
Since the operator ILL + n is non-negative and fi < we have that (IlL + n — /x) is invertible and 

Therefore the eigenvalue \i is simple. We now substitute v back into the eigenvalue equation for u and 
take a scalar product with We obtain that 

-{L + (UL + U - /x) _1 ILL + *, <f } + (L+tf , *) = fi{9, *) 

We now assume that there is another negative eigenvalue A. Then the corresponding equation also 
holds with A replacing fx. Subtracting we obtain 

-(l + (ul + u - / u)- 1 nL+^, *) + (L + (nx + n - A) _1 nL+*, *} = - a)(*, 

Using the resolvent identity 

(nx+n - - (nx+n - A)- 1 = (m - A)(nx + n - ^(nx+n - A) -1 

Thus, 

((nx+n - //) -1 (nx + n - A) _1 nL + \i/,nx + \i/) = -(*,*) 

Observe that (nX + n — //) _1 (nX + n — A) -1 is a positive operator since IlX + n is non-negative and 
A, < 0. Contradiction. □ 

14 Verification of the needed properties of L + and L_ for a particular 
class of nonlinearities 

The purpose of this section is to verify that the conditions imposed on L + and L_ in Definition 12.9 
hold for a particular class of nonlinearities. This class is chosen so as to allow for the nonlinear analysis 
of the previous sections to go through. 

(14.1) Ms 2 ) = s?- 1 /(g2) 



+ f(s 2 ) 

with a constant 9 > and the function / satisfying the conditions 

(14.2) C lS r+1 < \f(s 2 ) < C 2 s r +\ s 2 \f'(s 2 )\ < C 2 \f(s 2 )\, p e (1, 1 + -), r > -1 



n 



The existence of a ground state for the problem 



(14.3) _I A 0-/?(|0| 2 )0 + ^ = O 

for a/0 had been established by Berestycki and Lions under the following conditions on the function 
0: 
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1. > lim s ^ +00 /3(s)s~^ > +00 

2. There exists s > such that G(s ) = Jq° (3(s 2 )sds - ^s 2 > 
Moreover, in the case when the function (3(s) satisfies a stronger condition that 

(14.4) lim (3{s)s~~^ = 

a ground state can be constructed from a solution of the constrained minimization problem for the 
following functional: 

(14.5) J[u] = { [ \Vu\ 2 : W[u] = [ G(u) = l} 
If w is a minimum it solves the equation 

~Aw - \{f3{w 2 )w - yH)) = 

where the Lagrange multiplier A is determined from the condition that W[w] = 1. We can then find a 
ground state via rescaling 

(14.6) (f>(x) = w(X'^x) 

Observe that it is possible to choose w a positive spherically symmetric function. We now consider the 
case of the monomial subcritical nonlinearity (3(s) = with p < By the results of Coffman, 

McLeod-Serrin, and Kwong there exists a unique positive radial solution of the equation (14.3) for 
a^O. Let w denote the corresponding minimizer of the functional J. 

Definition 14.1. Given 7 > and w the minimizer of J corresponding to the unique ground state 4> 
define 

(7(7) = inf {9 : for any positive non-increasing radial function u with the 

(14.7) property that \\u — w\\ H i > 9 we have that J[u] > J[w] +7} 

We now make the following claim 

Lemma 14.2. Let f3{s 2 ) = s^ 1 be a monomial nonlinearity with p £ (1,^2). Then the function 
(7(7) — > as 7 — > 0. 

Proof. We argue by contradiction. Assume that there exists a sequence 7^ — ► 0, a positive constant 
a, and positive radial functions u lk such that \\u^ k — w\\ H i > 9 but J[w 7 J < J[w] + 7^. Then the 
sequence u lk ^ is minimizing for the functional J. This implies that 

||Vu 7 J L 2 - \\VW\\ L 2 

Using the constraint W^ 7 J = 1 it is not difficult to show that the sequence u yk is uniformly bounded 
in H 1 , see [BL]. Thus without loss of generality we assume that u~ jk — > u weakly in H 1 for some 
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\w\ p+1 



\w\ 2 



radial non-increasing function u. Therefore, u is another minimizer of the functional J and its rescaled 
version is a non-increasing radial solution of the equation (14.3). By the strong maximum principle 
it is positive 6 and therefore a ground state. Since the ground state is unique, after rescaling back we 
conclude that u = w. Therefore, we have constructed a sequence u lk with the properties that 

(14.8) u lk — > w weakly in H 1 , 

(14.9) Vu lk -► Vw in L 2 , 

(14.10) / 1^1^ = 1 + ^ 

(14.11) \\u lk -w\\ H i > e 

Since 2 < p + 1 < conditions (14.8) and (14.9) imply that 

/ Kr 1 - / 

Thus from (14.10) 

/ \ u Jk I 2 " 

and with the help of (14.8) and (14.9) we conclude that n 7fc — ► w in fl" 1 . This contradicts (14.11). □ 

14.1 Variational continuity 

We now consider the ground state problem 

for the nonlinearities 

(14.12) a^) = 

where 1 < p < 1 + ^ , the function / satisfies the following estimate 

Cis q < \f(s 2 )\ <C 2 si +1 
for some positive constants C\, C2 and q > — 1, and > is a positive number. Define 

(14.13) G e (r)= / /? e ( s 2 ) S( fe-— r 2 , 

io 4 

(14.14) W e [u}= f G e (u(x))dx 

6 The minimizer w cannot be identically zero since one can show that the minimum is attained on the function satisfying 
the constraint W[u] = 1, see [BL]. 
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Lemma 14.3. We have the following estimate 

(14.15) \Wg[u]-W [u]\<9^ f (|u| 2 + \u\ p+1 ) 

JK. n 

Proof. Estimate (14.15) immediately follows from the inequality 



(14-16) \G e (r) - G (r)\ < r p+1 ——, = fh*-*— 



3+1 



' + ri +1 e + t<? +1 

since the above expression can be bounded by 

min {9T p - q ,T p+1 } 

i i 
Thus using the first term for the values of r > 9p+i, and the second term when r < 0p+i we obtain 

(14.16) . □ 
We now consider the variational problem 

(14.17) Jg[u] = { [ \Vu\ 2 : Wg[u] = l} 

Proposition 14.4. Let <p be the ground state of the problem (gro). Then for any sufficiently small 
9 > there exists a positive constant 5' = S'(0) — > as 9 — > 0, and a ground state (fig of (grg ) such that 
\\cfig-cfi\\ m <5'. 

Proof. We start by choosing a sufficiently large constant M such that for all sufficiently small 9 any 
minimizer of Jg is contained in a ball B M / 2 of radius M/2 in the space H l . In particular, using (14.15) 
we will assume that for u € Sm 

(14.18) \Wg[u]-W [u]\<e^ 

We now observe the following trivial property of the constraint functionals for any 9 > and 

an arbitrary fj, ^ 

(14.19) Wg[u(x)] = l j n Wg[u{-)} 

We now fix a sufficiently small 9 > 0. Let w be the minimizer of the variational problem J = Jq 
corresponding to the unique ground state <fi. The function w satisfies the constraint Wo [if] = 1. 
Therefore, using the rescaling property (14.19) and (14.18) we can show that there exists \i = n(w) 
with the property that 

Wg[w{-)} = 1, 

(14.20) |/i-l|<0^5 
Moreover, 

(14.21) JeW-)} = V n - 2 J [w(x)} = JqH + 0(#FFi) 
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We now claim that there exists a small positive 5 = 5(9) — > as 9 — > 0, such that for any positive 
non-increasing radial function u satisfying the constraint Wg [u] = 1 and the property that 

(14.22) \\u-w(-)\\ H i > 5 

A 1 

we have 

(14.23) J e [u] > Jg[w(-)} + 9^ 

A 1 



Assume for the moment that the claim holds. Then (14.22) and (14.23) imply that Jg has a minimizer 
in th 
that 



in the S neighborhood of the function w(-). We denote this minimizer by wg. Then (14.20) implies 



\\wg - w\\ H l < \\wg - w(-)\\ H i + 1 1 10 - w(-)\\ H i <5+\\w- w(-)\\ H l 

Observe that \\w — w(-)\\jji — > as /x — > 1, which follows by the density argument and the fact that 
it is easily satisfied on functions of compact support . Define the function a w (e): 

x 

(14.24) a w (e) := sup \\w - w(—)\\ H i, a w (9) — > as 9 — > 

p-i A* 

|^-1|<6»pT9 

Therefore, 

(14.25) \\w e -w\\ H i<S + a w {9) 

The functions wq , to are the solutions of the Euler-Lagrange equations 

1 a 2 

(14.26) - - Aw 9 - \e{f3e{w e ?)w e - —w e ) = 0, 



(11.27) - ^Aw - \((3 (w) 2 )w - yw) = 



where the Lagrange multipliers Xg, X are determined from the conditions that W^u^] = Wo [to] = 1. 
We multiply the equations (14.26) and (14.27) by wg and w correspondingly, integrate by parts, and 
subtract one from another. Using the estimate 

/ \Pe(wg) 2 )w 2 g - (3o(wg) 2 )w 2 \ < 9^, 
JM. n 

which is essentially the same as the estimate (14.18), and the estimate (14.25) we obtain that 
(14.28) (A - A,) / (A)H> 2 - ^-w 2 ) = 0(5 2 ) + 0(9^) + a 2 w (9) 

Recall that (3q(w 2 ) = io p-1 . The condition that W[w] = 1 implies that 



/ 



p+ 1 4 1 1 ; 



7 In fact, the minimizer w is smooth and localized in space and thus one could even give the precise dependence on fj, 
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Thus, 

This allows us to conclude that 

(14.29) |A-A e | < <5 2 + fl^F 1 + a 2 w {9) 

Finally, recall that the ground states 4>g and (f> are obtained by the rescaling of the minimizers we and 

_ I 

(j)g(x) = Wg(X g 2 X), (f)(x)=w(X 2 X ) 

Thus 

_l _i 
\\<t>o ~ 4>\\hi ^ lkfl(A„ 2 x) - w(X ^x)\\ H i 

(— — — 1\ A# i 

x e + X e )\\ w o( x ) -w{x)\\ m + \\w(x) - w{{— )*x) \\ H i 

By (14.29) the constants A# are uniformly bounded in terms of the absolute constant A, which depends 
only on w. Moreover, Xe — > A as 6 — > 0. We appeal again to the H 1 modulus of continuity of the 
minimizer w and define the function 8 

(14.30) b w (6,8) := sup \\w(x) - w(fj,~^x)\\ H i 

\^-l\<C(S 2 +e^¥+i 1 +al(0)) 

for some positive constant C dependent only on the minimizer w. The function b w (9, S) — > as 0, S — > 0. 
Therefore, since we have already proved in (14.25) that w$ is close to to in H 1 , we obtain 

(14.31) U e - <f>\\ H i <5 + 9^ + a w {6) + b w (6, 5) 

Since by the claim S = 5(9) — > as 6* — > and the functions a w (6), b w (9, 5) also have this property we 
obtain the desired conclusion. 

It remains to prove the claim (14.22), (14.23). Let u be as in the claim, i.e., u <G Bm and Wg[u] = 1, 
and 

(14.32) \\u - w{-)\\ H i > 5. 
A* 

s One can show that a / the Lagrange multiplier A / 0. This follows from the following argument. By interpolation 
for p < 2+§ 

1 — n — 2 



J ^ +1 <||vhi2 2 V IM 



L2 



Thus for n > 2 the power p + 1 — n 2 ^ < 2 and using Cauchy-Schwarz, constraint = 1 and the assumption that 

a / 0, we can show that ||Vw|| i 2 > c for some positive constant c. Repeating argument determining the Lagrange 
multiplier we verify that A^0 
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for some 5 to be chosen below. Similar to (14.20) we can find a constant v = v{u) such that 

(14.33) WoK^)] = 1, JoK^)] = + O(0^?), 

(14.34) \v-l\<9^ 

Using (14.20), (14.32), (14.34), and definition (14.24) we infer that 

(14.35) \\u(-)-w\\ H i > ||u(-) -w(— )\\ H i - \\w(— ) -w(x)\\ H i >u^5-a w {6) >5-6^5-a w {0) 

We now use Lemma 14.2 for the variational problem J = Jq. This gives a function (7(7), with the 
property that (7(7) — > as 7 — > 0, such that for any radial non-increasing positive v with the property 
that \\v — w\\hi > a(-y) and WoM = 1 we have Jq[v] > Jq[w] + 7. We set 

7 = 56»i+i, 5 = a(j) +6^ +a w {9) 

It follows from Definition 14.1 of 17(7) and (14.35) that with these choices, function u(|) verifies the 
inequality 

JoH^)] > J [w]+5e** 

Finally, using (14.21) and (14.33) we obtain 

J e [u] > Je[w(-)] +36^ 

It remains to note that the constant 5 in has been chosen 

S = a(0^) + 6>?+i + a w {6) 
and by Lemma 14.2 and (14.24) goes to zero as 6 — > 0, as claimed. □ 
Recall definition of the operator L e + associated with the ground state <pg. 

(14-36) L% = ~A - p e ($) - 20 , ($)<t% + ^ 

Denote 

(14.37) V e = Pe(4>e) + ^M)4>e 

Using the definition of (3$ we compute V$ explicitly 

a2\ affj.2 



(1438) ^W + 2 ^ ^71 



Uniform bounds on ground states (f>g guaranteed by the Proposition 14.4 imply the following result. 
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Lemma 14.5. Let function f obey the assumptions that 

C lS i +1 < \f(s 2 )\ < C 2 s« + \ S 2 \f(s 2 )\ < Cif(s 2 ) 

for some positive constants C\, C 2 and q > — 1. Then for any p € (1, 1 + ^] there exists a r = r(p) in 
the interval r € [§, 00) such that 

(14.39) \\Vo-Vo\\ L r ^0, 0^0 

Proof. We have a pointwise bound 

1 V— 1 —1 2 2n 

In addition, since <fo — > in we have that <p e — > in the space Lp- 1 n L^^K^ 1 ) . Since 

' 2 2n i 

(n-2)(p-l)J 



we obtain the desired conclusion for any r in the interval r € (—§4, , _ 2 2 w The existence of the 



Lebesgue exponent r in the desired interval now follows from the restriction p£ (1,1 + -] on the range 
of the exponent p. □ 

Corollary 14.6. The operators 

(4-L+X-A + l)- 1 : L 2 ^L 2 , 
(-A + l)- 1 (L e + -L + ) : L 2 — ► L 2 

w/if/i f/ie norm converging to as # — > 0. 

Proof. The difference — L + = V# — Vo- The result now follows from Lemma 14.5, Sobolev embed- 
dings, and Holder inequality. □ 

14.2 Stability of the convexity condition 

The goal of this section is to prove that the ground states 4>e of the problem 

1 a 2 

-A0 + /3 e (0 2 ) = T 

with the nonlinearity 
(14.40) /W 2 ) = M p -V: 



e + f(<py 



C lS i +1 < \f(s 2 )\ < C 2 s* + \ s 2 \f'(s 2 )\ < C 2 \f(s 2 )\, p € (1, -), q > -1 



n 



verify the monotonicity condition (1.17), {L e + 1 (/>o, 4>q) > 0. 
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Theorem 14.7. Let <f>g be ground states constructed in Proposition 14-4 with the nonlinearities (3q 
satisfying (14.40). Then then for a given nonlinear eigenvalue a / and all sufficiently small 9 the 
ground states 4>q verify the monotonicity condition 

(14.41) (L d + ~ 1 4> e ,4>e)>0 

Proof. Condition (14.41) is meaningful provided that 4>q is orthogonal to the kernel of L d + . We start 
by examining the spectrum of the operator L + . The operator L + has a unique negative eigenvalue, 
the zero eigenvalue has multiplicity n and the corresponding eigenspace is spanned by the function 
7^4>, [We2]. The rest of the spectrum is contained in the set [\, oo). Therefore, in the case a / the 
spectrum S(L + ) of L + has an isolated discrete component (in fact two components). We can construct 
an eigenspace projector Pq of an isolated component of the discrete spectrum 

(14.42) ft-JLjjW-,)-* 

with an arbitrary curve 7 encircling the desired spectral set and such that 7 PI S(L + ) = 0. Consider 
now the resolvent of at z such that dist(z, X(L + )) > C for some sufficiently small constant C, 
which only depends on L + . We have 

(14.43) {L% - z)- 1 = (L + - z)- 1 - (4 - z)~\L% - L + )(L + - z)' 1 
It is not difficult to show that for such z 

\\{L + -z)- l f\\ m <\\f\\ L , 
Therefore, using Corollary 14.6 we can conclude from (14.43) that for all sufficiently small th > 

||(Lt-^)- 1 ||<2||(L + -z)- 1 || 

and thus z ^ S(L^_). Moreover, 

(14.44) || (L e + - z)" 1 - (L+ - z)' 1 1| < c(0) 

for any z : dist(z, £(L + )) > C. By Corollary 14.6 the constant c(9) — > as — > 0. Therefore, for the 
same path 7 as in (14.42) we can define 

(14.45) fl-JLjf <L» -,)-* 

Moreover, for all sufficiently small 9 > the rank of Pq remains constant. Thus, for any sufficiently 
small 6 the operator L e + has a unique simple negative eigenvalue and a zero eigenspace of dimension 
n. Since we know that the functions ^r-<f>9 are contained in that subspace, they, in fact, span it. 
Therefore, 4>e is orthogonal to the kernel of L e + and the expression (14.41) is well defined. 

For any sufficiently small 9 > we set Qg to be a projection on the orthogonal complement of the 
null eigenspace of L 6 + . Let A U e S(L^). Define the operators 

(14.46) Kg(X) := Q O (L + - \)~ 1 Qo - (L + - A)~ 1 Q 
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It follows from (14.44) and the properties of the spectrum of L + that for all small 6 > and all A such 
that |A| < C 

(14.47, ^-^Q4< Aist{x J Llmo)) <C- 

for some universal constant C", determined by the operator L + . Also note that 

(14.48) \\Qe-Qo\\<c(9) 
This is a consequence of (14.44) and the definition 

with a short path 7 around the origin. Using the resolvent identity 

{L% - z)~ l = (L + - z)' 1 + (L + - z)-\L e + - L + )(L e + - z)~ l 

we obtain that for any A UeS(L^_) 

K e {\) =Qo(L + - X)~ l Q e - (L + - A)' 1 Q + Qo(L+ - \T\L e + - L + ){L e + - X)' l Q e 
=(L+ - X^QoiQe - Qo) + (L+ - A)- 1 Q (4 - L+){L% - \)- l Q e 

Using Corollary 14.6, (14.47), and (14.48) we infer that for any A < c and A UeT,(L 9 + ) 

(14.49) \\K e (X)\\ < c{9){\ + \\{L + - A)- 1 Q (-A + 1)||) < c'(9) 

uniformly in A. The last inequality follows since the operator norm of (L+ — A) _1 Qo(— A + 1) is 
bounded by a universal constant dependent on L + only. This can be seen as follows. Since Vq is a 
smooth potential and (L + — X)~ 1 Qo is bounded on I? we can replace the operator (— A + l) by (L + — A) 
and the result follows immediately. 

We now test the operator Kg{X) on the ground state 4>g. Using that Qo<po = 4>q we obtain 

K e (X)(j> e = Q (L e + - A)- Vfl - (L+ - A)- V + (L+ - A)- 1 Q (</> - <j>e) 
Coupling the above identity with (p. 

{fa, (L e + - A)" Vfl) " (<t>, {L+ - A)"V> = K e {X)4> e ) + ((4>e - <h), {L% " A)~V*> 

+ (0, (L + -A)- 1 Q o (0-«)=O(c(0)) 

where we have used that Qe<Pe = 4>0i the bound (14.47), and the estimate \\<pQ — </>||#i < c(8), which 
follows from Proposition 14.4. The above holds uniformly for all |A| < c and A U#£(L+). Passing to 
the limit A — > 0, say from the upper half-plane, we obtain that for all sufficiently small 9 > 

(4> e , (L 6 ,)- 1 ^) = (0, LlV> + O(c(0)) < 



u -i i - i 1 \1 * 1 rvi-i \i 1 1 1 * 1 u I ■ 1 I /1 1 n / 



The last inequality follows since by the assumption is a stable ground state, i.e., {<fi, L, 1 ^) < 0. □ 
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